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ABSTRACT 


A physical understanding and consequent mathematical 
modeling of RF energy in naval indoor environments is of vital 
importance to the usability and effectiveness of communication 
systems used by Navy. Over the last few years, there is a 
growing interest in placing Wireless Local Area Networks (WLANs) 
in ships and submarines. Especially large ships yet to be 
constructed, are designed with increased electronic systems but 
limited personnel. Reliable electronic systems will be crucial 
for efficient ship operation and survivability. 

This thesis investigates the feasibility of deploying a 
physical model called Numerical Electromagnetic Code-Basic 
Scattering Code (NEC-BSC) to simulate confined naval compartments 
in the 2.4 GHz Industrial Scientific Medical (ISM) band. More 
specifically, the coverage area, the number and positions of 
transmitters and observation points and the statistics of Radio 
Frequency (RF) signal distribution were described. The area 
specifically targeted for this research was a typical two-story 
missile room. Additionally, some important conclusions regarding 
the validity of NEC-BSC for indoor applications are presented and 
some recommendations for future research are provided. 
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EXECUTIVE SUMMARY 


Indoor wireless systems are already used in a huge 
variety of office, residential and even military 
environments. Two types of models, the empirical and the 
physical ones, are powerful tools, capable of simulating a 
given confined area and collecting significant databases. 

The Numerical Electromagnetic Code-Basic Scattering 
Code (NEC-BSC), developed at Ohio State University under 
government contract, is a well-known physical model in 
electromagnetic field theory. Over the last ten years it 
has been of great use in analyzing scattering and radiation 
problems and visualizing the flow of energy from antennas. 
Given that the code's applicability has recently been 
proved reliable for complex buildings [37] , we tested its 
performance to a typical missile room of a military ship in 
the 2.4 GHz ISM band. 

At first, the case of the electric field returned by a 
smooth perfectly conducting convex surface was 
investigated. Then, the results•were compared to the ones 
obtained by the reliable software called "PATCH", within 
the framework and limitations of the Method of Moments 
(MOM) . The comparison offered encouraging results. In 
addition, other simple geometries were simulated in order 
to gain familiarity with the code and its limitations. The 


xi 



goals, however, were to construct a typical missile room 
with material (metallic) uniformity, examine the electric 
field strength coverage maps and perform a statistical 
analysis of the electric field strength variation. 

To this end, the antenna was placed in four different 
positions and data were collected. Then the excess "room 
gain" over free space of the same volume was calculated 
within certain distances from the antenna, called "distance 
groups." Next, the probability density function (pdf) and 
cumulative distribution function (cdf) of the variation of 
the "room gain" variable was introduced. It appeared that 
the key parameter in defining the mean value of the excess 
gain was the relative number of observation points that 
were shadowed by the obstacles and thus experience non- 
line-of-sight propagation attenuation. In any case, 
however, the distribution followed the lognormal one and 
its pdf and cdf obeyed the Gaussian ones. Finally, the 
mean value of the "room gain" ranged from -0.1 dB to + 2.91 
dB. 
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I. 


INTRODUCTION 


A. WIRED VS WIRELESS COMMUNICATION 

The last century has seen the rise of telephony, as 
well as the successful deployment of wireless 
communications. • The introduction of these inventions 
underscored, without doubt, a fundamental change in the way 
people communicated and was a critical step toward 
shattering the barriers of time and space. To date, 
however, the hindrance of location has not been conquered, 
since most individuals are still confined to wired 
equipment for communications [1]. 

Rapid changes in word communications have been 

prompted by demand for more reliable and convenient ways to 

communicate. Especially during the past decade, wireless 

services proved to be cost effective solutions for 

providing data due to their inherent advantages over the 

conventional, wired services. These advantages include 

simplicity and ease of use, integrity and reliability, 

compatibility with existing networks without the need for 

expensive and time-consuming rewiring. In addition, there 

is more immunity to natural and man made disasters, 

flexibility of shifting terminals and equipment around, 

relatively easy maintenance, providing services for a 
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larger number of users, roaming and, finally, elimination 
of wiring around and inside buildings. The introduction of 
successful wireless systems is comparable in significance 
with that of any major invention made during the 20th 
century. The deployment of wireless systems is increasing, 
while the wired communications is holding steady if not 
decreasing. 

All the aforementioned advantages of using wireless 
communication mainly in small or large offices, public 
services buildings and even military environments, have 
built a potential market of a very large number of users. 
In order to meet these demands, new systems have been and 
are being designed to provide services of the best possible 
quality in indoor and outdoor environment [2]. 

B. WIRELESS PROPAGATION BACKGROUND 

Radio communication, to date, is enjoying the fastest 
growth period in history, since enabling technologies 
permit wide spread deployment. During the past decade, 
radio communications business has grown by orders of 
magnitude, fueled by large-scale circuit integration, 
digital and radio frequency (RF) circuit fabrication 
improvements and other miniaturization technologies which 
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make radio equipment smaller, less heavy and more 
affordable and reliable [3]. 

The Ultra High Frequency (UHF) electromagnetic 
spectrum ranges from 300 MHz to 3 GHz. There is much more 
frequency space available for communications in this band 
than in the entire spectrum below. Moreover, high 
directivity obtainable at microwave frequencies makes reuse 
of these frequencies possible.in the same area. 

The most basic model of radio wave propagation 
involves the ideal "free space" model which treats the 
region between the transmit and receive antennas as being 
free of all objects that might absorb, refract, reflect, 
diffract or scatter radio frequency energy. It thus 
follows that, within this region, the medium behaves as a 
perfectly homogeneous and nonabsorbent medium. 
Furthermore, earth is treated as being infinitely far away 
from the propagating signal (or, equivalently, as having a 
reflection coefficient that is negligible). 

In free space the spreading of RF energy creates 
apparent attenuation between the transmitter and receiver 
according to an inverse-square law. In typical 
broadcasting applications radio waves emanate from a point 
source of radio energy, traveling in all directions along 
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straight lines, filling the entire spherical volume of 
space with radio energy that varies in strength with a 
1/(range) 2 rule (or 20 dB per decade increase in range). 
The received power, expressed in terms of transmitted 
power, is attenuated by a factor called (free space) path 
loss. 

Calculating free space loss is quite simple. Consider a 
transmitter with power Pt coupled to an antenna, which 
radiates equally in all directions (an isotropic antenna). 
At a distance d from the transmitter, the radiated power is 
distributed uniformly over an area of 4d 2 (i.e.’the surface 
area of a sphere of radius d), so that the power flux 
density is: 

s = —, Watts/m 2 . (1.1) 
And 2 

The transmission loss then depends on how much of this 
power is captured by the receiving antenna. If the capture 
area, or effective aperture of this antenna is Ae then the 
power, which can be delivered to the receiver (assuming no 
mismatch or feed line losses), is simply: 
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P r = sA e Watts. 


(1.2) 


For the hypothetical isotropic receiving antenna and a 
wavelength X we have: 


A 


e 



(1.3) 


Combining all equations, we obtain: 


P R =P, 


And 


Watts. 


(1.4) 


The free space path loss between isotropic antennas is 
then Pt/Pr. For this case of idealized propagation, received 
signal power is very predictable. 

In reality, however, propagation rarely follows the 
free space model. Three major mechanisms, which are going 
to be examined in the next pages, dominate in the radio 
propagation: (1) reflection, (2) diffraction and (3) 
scattering. All three of these phenomena result in radio 
signal distortions, additional losses, and cause signals to 
fade. 
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C. THE PHYSICS OF PROPAGATION 

The phenomena, which govern radio propagation, can be 
very complex and diverse. However, we can categorize them 
according to the three major mechanisms of object 
interaction (reflection, diffraction and scattering). 

Reflection occurs when a propagating electromagnetic 
wave impinges upon an object, which has very large 
dimensions compared to the wavelength of the propagating 
wave. Reflections occur from the surface of the earth and 
from buildings and walls, and reflected waves can interfere 
constructively or destructively at the receiver [3]. 

Diffraction occurs when the radio path between the 
transmitter and receiver is obstructed by an impenetrable 
surface. The secondary waves resulting from the 
obstructing surface are formed throughout the space and 
even behind the obstacle, giving rise to an apparent 
bending of waves around the obstacle, enabling reception 
even when a line-of-sight (LOS) path does not exist between 
transmitter and receiver. Diffraction is assumed to be a 
strictly local phenomenon; this means that the diffracted 
field depends only on the power of the incident field at 
the exact diffraction point, and on the local geometry of 
the diffracting object [3]. 
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At high frequencies, diffraction, like reflection, 
depends on the geometry of the object, as well as the 
amplitude, phase, and polarization of the incident wave at 
the point of diffraction. Because of diffraction RF energy 
can be received, eventhough the source is "shadowed" by an 
obstacle [3]. 

Scattering occurs when the medium through which the 
wave travels contains objects ("obstacles") with dimensions 
that are small compared to the wavelength, and where the 
number of obstacles per unit volume is large. Scattered 
waves are produced by rough surfaces, atmospherics (rain 
and snow), or by other irregularities in the channel. 
Foliage, street signs, and lampposts induce scattering in a 
mobile communications system. Figure 1 shows the different 
propagation mechanisms in the presence of an arbitrary, 
impenetrable obstacle. 
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multiple signal propagation paths exist, caused by whatever 
phenomenon, the actual received signal level is a vector 
sum of the all signals incident from all directions. Some 
signals will aid the direct path (assuming that the direct 
path exists), while other signals will subtract (or tend to 
vector cancel) from the direct signal path. The total 
composite phenomenon is called multipath. 

As a mobile receiver moves through a specific indoor 
space, the three phenomena affect the instantaneous 
received signal in different ways. For a clear LOS path to 
the transmitter, diffraction and scattering will not 
dominate. On the other hand, in an environment where there 
are many obstacles of different sizes and constituent 
materials and the direct channel is blocked between the 
receiver and the transmitter, diffraction as well as 
scattering are most likely to dominate [5]. 

Two kinds of multipath exist: (1) specular multipath 
arising from discrete, coherent reflections from smooth 
metal surfaces, and (2) diffuse multipath arising from 
diffuse scatterers and sources of diffraction (the visible 
glint of sunlight off a choppy sea is such an example) . 
Both forms of multipath are undesired for radio 
communications. Diffuse multipath provides a sort of 
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background "noise" level of interference, while specular 
multipath can actually cause complete signal outages and 
radio "dead spots" within a building. 

Another phenomenon, called scintillation, describes 
the rapid fluctuation of the phase and intensity of a radio 
signal that has passed through the earth's ionosphere, 
typically on a satellite-to-ground propagation channel. 
This phenomenon is similar to the twinkling of stars in the 
night sky, except in this case the fluctuations are caused 
by small-scale variations (irregularities) in the 
ionospheric plasma density along the propagation path 
followed by the signal. 

D. FUNDAMENTALS OF ANTENNAS 

1. Radiation Pattern 

The antenna's radiation pattern is the graphical 
representation of the radiation characteristics of the 
antenna as a function of space coordinates. The most usual 
representation is the three-dimensional spatial 
distribution of radiated energy along an observer's 
position at a constant distance from the antenna. The 
graph of the spatial variation of the electric field along 
that path is called a field pattern. 
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An isotropic antenna is defined as a hypothetical 
lossless antenna that radiates or receives equally in all 
directions. Clearly, an isotropic antenna is a fictitious 
entity, since even the simplest antenna has some degree of 
directivity. However, an isotropic radiator represents a 
convenient reference for comparing the directional 
properties of physical antennas. On the other hand, a 
directional antenna radiates or receives electromagnetic 
energy more effectively in some directions than in others. 
As an example, the radiation pattern of a half-wave linear 
dipole is shown in Figure 2. 



Figure 2. Power Pattern of a Half-Wave Dipole {from [6]). 
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2. 


Near Field and Far Field Regions 


The antenna's electromagnetic (EM) field can be 
classified according to the relative (in terms of operating 
wavelength) distance to the antenna as: 

• The reactive near field region adjacent to the 
antenna. In this case, the field patterns change 
rapidly and reactive energy is stored around the 
antenna 

• The radiating near field region (for large 
antennas known as the Fresnel region). Again, 
the field patterns change rapidly, but radiating 
energy oscillates towards the antenna. 

• The furthermost, the far field region, commonly 
referred to as the Fraunhofer region, has no 
reactive energy and is represented by a spherical 
wave (for any realistic antenna of finite size) 
propagating away from the antenna. 

Certain behavior characteristics of EM fields dominate 
at one range from the radiating antenna, while others, 
completely different, can dominate at another distance. In 
order to categorize characteristics of electromagnetic 
fields as a function of distance from the radiating source, 
we define the following boundaries: the "Near-Field 
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boundary," the "Far-Field boundary" and the "Transition 
Zone boundary," which has a combination of the 
characteristics found in both fields. The regional 
boundaries are usually categorized as a function of the 
operating wavelength (X ) . 

We must state, however, that these regions categorize 
behaviors, which may vary even within each region. In 
addition, the boundaries for these regions are approximate 
"rules of thumb" (more precise boundaries can be defined 
based primarily on antenna type and antenna size, and even 
then the experts differ). Figure 3 shows these regions and 
boundaries. 
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Figure 3. Antenna Field Regions for Typical Antennas (from 
Ref. [ 7 ] ) . 

When placing an antenna, it is of vital importance to 
remove all objects from at least the near field, since 
their presence can alter the design radiation and impedance 
characteristics unpredictably. 

E. OBJECTIVE AND APPROACH 

The objective of this thesis is to evaluate the 

Numerical Electromagnetic Code - Basic Scattering Code 

(NEC-BSC), for shipboard (indoors) wireless propagation 

prediction. Using the code, we will predict the signal 

strength distribution for wireless communication in the 2.4 
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GHz- Instrumentation-Scientific-Medical (ISM) band for a 

given antenna location in a shipboard compartment. 
Simulation results will be presented graphically and the 
signal strength statistics will be extracted from the 
results. 

F. THESIS OUTLINE 

This thesis consists of six chapters and is organized 
as follows. Chapter I compares the wired and wireless 
communication and provides the introduction of the RF 
propagation. Chapter II describes the theoretical 
principles of indoor propagation. Chapter III presents the 
NEC-BSC software overview. Chapter IV describes some basic 

examples and the associated results obtained by using NEC- 
BSC code. Chapter V provides the development of a 
realistic model of a shipboard compartment, the simulation 
results and the statistics for signal propagation. 
Finally, chapter VI concludes the thesis and present 
recommendations for future work. 
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II. INDOOR PROPAGATION MODELS AND THESIS OBJECTIVES 


A. COMPARISON BETWEEN INDOOR AND OUTDOOR MODELS 

Outdoor environments are characterized by significant 
signal strength fluctuations caused by reflection, 
diffraction, and scattering. These signal strength 
variations may have a range of 30 to 40 dB at Ultra High 
Frequencies (UHF) and account for some of the difficulties 
that a designer of reliable radio communications systems 
has to cope with. For example, if we shift the receiver by 
a distance of a half-wavelength in a random direction, the 
received signal magnitude may change drastically. The 
basic signal attenuation with range noticed in the real 
world describes the so-called "large scale" effects, while 
the signal strength fluctuations with motion are reported 
"small scale” effects [8] . 

Outdoor propagation models presented by Ibrahim [Ref. 
10], Okumura [11] and Hata [12], have been developed for 
large cells, whereas for indoor communications the suitable 
cell diameter is substantially smaller and such cells are 
usually termed as nanocells and picocells. Hence, the 
outdoor models are inadequate for accurate indoor 
propagation prediction. 
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Indoor propagation, although dominated by the same 
mechanisms, can not be treated as an extension of outdoor 
propagation. For example, signal levels depend on whether 
interior doors and windows are open or closed. Also, the 
location of antennas affects the large-scale phenomena, 
since an antenna mounted at a floor level in a partitioned 
room receives different signals than that same antenna 
mounted above the floor. Moreover, in a room crowded with 
obstacles, all propagation mechanisms become more and more 
complex and the .smaller propagation range makes it even 
more difficult to ensure far-field radiation for all 
plausible types of antennas and receiver positions. 

B. MODELS FOR INDOOR PROPAGATION 

To make things more challenging, it is impossible to 
design an "RF friendly" building, free of multipath 
reflections from numerous smooth surfaces, diffraction 
around sharp corners or scattering from walls and ceilings. 
In addition, the RF signal is subject to temporal 
variations due to the motion of people, doors, etc. [13], 
[14], [15], [16]. 

Radio wave propagation inside metal-walled 
compartments such as a missile room of a military ship is 
dominated by strong reflections off the compartment walls. 
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strength at 2.4 GHz can vary 


In such an environment, signal 
dramatically over a distance of a few centimeters. Despite 
these problems, we are still capable of predicting radio 
propagation in a typical room or a shipboard compartment. 
To this effect a number of models have been developed and 
semi-empirical formulas for indoor signal propagation 
prediction developed from propagation studies. Figure 4 
illustrates the results of a study of a signal at 1.8 GHz, 
penetrating a 'window? of a typical room. 



Figure 4. The Penetration of a Mobile Phone Signal at 1.8 
GHz Through a Window in a Reinforced Concrete Wall is 
Studied by Means of a Two-Dimensional Model. Right: 

Geometry. Left: Field Levels in dB Due to the Line Source 
in the Left Upper Corner. (From [17]). 

Generally, models can be classified as either 
statistical, based on empirical data, or deterministic, 
valid for the soecific environment 


under consideration. 









Some models do not belong in either category, since they 
are a combination of those two. These are referred to as 
hybrid models. 

Statistical models [18] and [3], are based on fitting 
curves or analytical expressions to observed empirical 
data, such that the propagation mechanisms are accounted 
for in an average sense. Since reality seldom corresponds 
to relatively simple mathematical models, a certain range 
of error cannot be avoided. The accuracy of these models 
depends on the accuracy of the measurements taken, as well 
as che similarities of the room to be examined with the 
(usually simpler) room in which the results were collected. 
However, the statistical models are simple to implement and 
computationally efficient. 

With Deterministic models, ('also called site specific 
(SISP) models [19], [20], [21], [22], [23], [24], [25] 
specific environments are characterized individually 
(rooms, buildings, tunnels, shipboard compartments etc.). 
In this way the model accuracy is improved, as long as 
there is detailed database available for the geometry, 
materials, location/motion of people etc for the indoor 
environment under consideration. The obvious and very 
important drawback in this case is that even if the 
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geometry/media database is available, the algorithms for 
indoor propagation calculation are typically very complex 
and, depending on the accuracy desired, model calculations 
can be extremely time consuming. 

In general, various "ray tracing" algorithms are 

employed for indoor propagation prediction. Such 

algorithms typically implement an exhaustive search of a 
"ray tree" as follows: first and foremost, a ray tracing 
code detects if there is a LOS path between the receiver 
and the transmitter and determines the RF energy received 

via the direct path. Then, depending on the source ray 

algorithm, the program detects if an object intersection 
occurs. If such is the case, code distinguishes the 
phenomenon as scattering or diffraction or reflection and 
computes the received signal, otherwise process stops and a 
new ray is "launched." Typically, as the number of rays 

encountered in the problem increases linearly, the number 
of required computation increases exponentially. 
Therefore, for given computational resources, there is a 
limit on the relative size of the indoor geometry 'that a 
particular code can be applied to. 

Hybrid models [19], [26] can be a successful 

combination of the statistical and site-specific models by 
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trading accuracy for computational efficiency. A survey of 
field strength prediction models is presented in [27]. 

NEC-BSC utilizes ray-tracing technique based on the 
Uniform Theory of Diffraction (UTD) and the geometry and 
materials data for the specific environment to be analyzed. 
Therefore it can be classified as a physical model. 

C. POTENTIAL BENEFITS FOR THE NAVY 

As mentioned previously, the majority of materials 
used in a shipboard compartment are metals (good 
conductors) that can be approximated as Perfect Electric 
Conductors (PECs). However, the geometry of a shipboard 
compartment (on the scale of the operating wavelength) can 
be very complex and the behavior of the transmitted signals 
can be difficult to predict. Furthermore, very limited and 
often classified research and documentation is available 
for shipboard wireless propagation. 

In 1998, a bistatic measurement system was developed 
for the 0.8 to 3 GHz frequency range [28]. Experiments 
took place onboard the former USS SHADWELL for a variety of 
types of confined rooms, locations of transmitters and 
receivers, as well as open or .closed door configurations. 
The results led to the conclusion that the expected 
elimination of RF energy due to the ship's metallic 
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structure (bulkheads, tanks, overheads, watertight doors, 
etc.) did not apply. The major phenomenon observed, was an 
approximate 10-dB signal attenuation for each closed 
passageway door and/or bulkhead [28]. 

Let us now examine the motives that urged the Navy to 
sponsor these and other similar measurements. Over the last 
few years, there is a growing interest in placing Wireless 
Local Area Networks (WLANs) in ships and submarines. Future 
generations of large ships are being designed with more and 
more electronic systems but with fewer personnel. 
Currently, it is generally agreed that WLANs are more 
efficient than conventional LANs. First, WLANs can provide 
services almost everywhere on the ship, without the 
physical limitations that cables have. WLANs also mitigate 
problems due to physical damage to wires or fiber optic 
cables. Finally, a sophisticated WLAN can also offer 
installation flexibility, reduced cost and ability to 
transfer vast amounts of sensor data related to condition, 
operation and damage control of the majority of ship's 
departments. 

D. OPERATION, IMPLEMENTATION AND SAFETY ISSUES 

So far, the most common and convenient application has 
been the hand-held walkie-talkies that crewmembers use. 
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During state of readiness, communication can be established 
from various locations even if all intermediate watertight 
doors are closed. However, this application could lead to 
mild to severe Electromagnetic Interference (EMI) problems 
[28] . 

For all potential WLAN applications, EMI is the 
system "bottleneck." Therefore, WLANs have to be designed 
using the best combination of transmitter power, receiver 
performance and operating frequency. In 1985, the Federal 
Communication Commission (FCC) authorized license-free 
spread spectrum RF communications for low power (less than 
1 watt) transmitters in the 910 MHz, 2.44 GHz, and 5.8 GHz 
ISM bands [29] . This enabled a rapid growth of spread 
spectrum WLAN industry, primarily for commercial 
applications. Lately, the technology is being applied to 
military use as well. Finally, since the transmitter power 
levels of WLAN systems are typically less than 1 Watt and 
radiated power density is attenuated at least 20 dB per 
decade increase in range, crewmembers present in a WLAN 
coverage area are not at risk to high RF exposure. 
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Ill. NEC-BSC SOFTWARE OVERVIEW 


A. INTRODUCTION TO NEC-BSC 

NEC-BCS, developed at Ohio State University under 
government contract, has been used over the last ten years 
for solving EM scattering and radiation problems [34]. 
More specifically, this user-oriented code is capable of 
analyzing near and far field radiation and scattering in an 
arbitrarily shaped room/compartment at UHF (300 MHz to 3 
GHz) and above. Rooms and compartments can be modeled by 
combining a number of basic structural shapes: multi-sided 
flat plates and finite elliptic curved surfaces like 
cylinders, cone frustums, and finite composite ellipsoids. 

Over the past years, NEC-BSC has been successfully 
used to predict scattering from superstructures of ships, 
wings of aircraft, large tanks and indoor compartments of 
space stations, simulated using perfectly conducting 
materials. The newest version of NEC-BSC (4.2) 

incorporates a thin dielectric slab option, hence the code 
can be used to model propagation within buildings with non- 
metallic walls [38]. 

B. NEC-BSC MODELS AND LIMITATIONS 

NEC-BSC utilizes advanced formulations based on the 
Geometric Theory of Diffraction (GTD) and Uniform Theory of 



Diffraction (UTD) [30], [31] to solve problems involving 
electromagnetic coupling, scattering, diffraction and 
computation of near- and far-field antenna patterns. 

Geometrical Theory of Diffraction is a high-frequency 
technique suitable for the analysis of antennas that are 
large in terms of wavelength (□) and located in a complex 
environment. Keller [32], [33] originally formulated GTD 
as an extension of Geometrical Optics (GO) . GTD proved to 
be very reliable in calculating the non-zero fields m a 
shadow region. This was accomplished by introducing a set 
of diffracted rays analogous to the reflected and 
transmitted rays of GO. 

In NEC-BSC, ray-tracing algorithms define the 
parameters of the field incident onto various obstacles. 
The diffracted rays are then defined by UTD for all 
incident rays, and added to the geometrical optics terms 
for the observation points of interest. Since most of the 
scattered and diffracted rays will encounter obstacles 
themselves, a rapidly increasing number of higher order 
scattered and diffracted rays will be created. In order to 
limit the number of computations for a solution with 
acceptable accuracy, a sophisticated code is required to 
"trace out" the various possible combinations that will not 
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contribute significantly to the total field at the selected 
observation points. In this manner only the components 
with the strongest contributions need to be calculated, 
thus increasing the code's computational efficiency. 

There are, however, some limitations of this user- 
oriented code, which basically occur due to the very nature 
of the algorithms employed [34] . Since all diffracted 
components are calculated using UTD, each plate should have 
edges of at least a wavelength long. In terms of curved 
surfaces their major and minor radii and length should be 
at least a wavelength in extent. Furthermore, for 
dielectric (non-metallic) obstacles, the source has to be 
located at least a wavelength away from the surface and the 
incident field angle of incidence should not be too close 
to the grazing angle. Also, the dielectric layer must be 
thin, at most half of the operating wavelength. For 
example, at 2.4 GHz (12.5cm wavelength) dielectric 
thickness must be less than 6.25 cm. Moreover, antenna 
elements have to be at least half a wavelength from all 
edges. In any case, the antenna can not be mounted on 
curved surfaces. 

Limitations on the NEC-BSC accuracy also stem from the 
high order term elimination. Eliminating higher order 
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diffracted rays improves computational efficiency, but it 
also makes the code unsuitable for certain types of EM 
calculations. For example, higher order diffracted rays 
are important for accurately predicting the radar cross 
section of a complex structure. Furthermore, flat plates 
and simple curved surface models are not always adequate 
for modeling complex environments found within buildings or 
shipboard compartments. Due to all these limitations, NEC- 
BSC program proved to give accurate results for 
approximately the first 25 dB of the results dynamic range. 

As a whole, the present NEC-BSC execute module will 
accommodate models with the following parameter limits: 


Pattern points 

— 

1801 

Defined pattern points 

= 

181 

Defined frequencies 

~ 

181 

Plates 

= 

200 

Plate edges 

= 

30 

Curved entities 

= 

200 

Curved sections 

= 

10 

Material types 

= 

25 

Material layers 

= 

20 

Sources 

= 

300 
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• Receivers 


300 


Antenna data points 

= 91 

Wires 

= 200 


C. MODELING THE STRUCTURE 

Most of the structures of practical interest can be 
successfully approximated by a combination of multi-sided 
flat plates and some finite elliptic surfaces. These can 
consist of either good conductors (even PECs) or 
dielectrics. A group of flat plates, each having a certain 
number of edges, can model most of the objects with slowly 
varying radii of curvature. In either case there is a 
certain procedure that one has to follow in order to 
construct an arbitrary object. 

If a plate has to be modeled, one first has to specify 
the corners in a three dimensional Cartesian coordinate 
system. The only condition is that all these corners 
should form a flat plate. If the points are difficult to 
specify, the user can simply construct the plate in two 
dimensions and then rotate or translate it into its desired 
location in three dimensions using the appropriate commands 
provided by the software. 

A plate can be attached to another plate in NEC-BSC in 
two different ways. First, two plates can share a pair of 
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corners, forming a wedge. The angle of the wedge is 
defined by the angle of the two plates, based on the 

direction of antenna's illumination. Alternately, one 
plate can be mounted slightly within the volume of the 
other. By piecing together a number of such isolated 
plates, a larger, more complex general structure can be 
approximated with sufficient accuracy. 

NEC-BSC also has the capability to implement 
electrically thin dielectric layers in order to model 
structures like non-metallic doors, walls, glass windows, 
furniture, etc. For flat dielectric plates, the user has 
to specify the thickness (not exceeding one half of the 

operating wavelength), relative permitivity, relative 
permeability, and the loss tangents for the material. 

If a surface with a non-trivial curvature has to be 
modeled, one can use elliptic cylinders, cone frustums and 
composite ellipsoids. In NEC-BSC all curved surfaces, are 
defined by their center coordinates and their major and 
minor axes. In case of an elliptic cylinder, the position 
of the first and the last cap is given with respect to the 

cylinder axis. The cone frustum is characterized by the 

position of its rim junctions. Finally, the ellipsoid is 
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constructed by setting different values of radii curvature 


along its maximum and minimum z-axis positions. 

D. MAJOR COMPONENTS OF NEC-BSC MAIN PROGRAM 

The NEC-BSC code is constructed in a systematic 
manner, meaning that all procedures are categorized and 
executed in modular sections. Figure 5 shows the major 
components of the NEC-BSC program. 


[initiali ze Data ] 



Figure 5. Block Diagram of NEC-BSC Main Program (from 
[34] ) . 
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In the input part the problem geometry is specified. 
The second section prepares all necessary parameters for 
later use. In the following section the fields resulting 
from individual propagation mechanisms are calculated and 
added to form the solution. Finally, all the outputs are 
provided to the user, either in the form of numerical 
results or as ray path visualization. 

E. TABLE OF COMMANDS 

A listing of NEC-BSC commands is provided in the 
preliminary draft of the NEC-BSC user's manual. (The 
listing of commands is still incomplete at the time of this 
writing) . The commands can be grouped into six broad 

categories: the geometry control commands, the structure 

definition commands, the antenna definition commands, the 

pattern control commands, the output control commands, and 

the program control commands. In creating a NEC-BSC input 
file the order of groups should be as given above, but 
there is no required order for the commands within the same 
group. The default unit is meters for all lengths and the 
default frequency is 299.9 MHz. All angles are specified 
in degrees. The following examples illustrate most of 
these points and validate a sample NEC-BSC solution. 
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IV. 


DEVELOPMENT OF BASIC EXAMPLES 


In this chapter, NEC-BSC results for several sample 
problems are presented. All pertinent Matlab, NEC-BSC, and 
Mathcad files are included in Appendices A through C, 
respectively. 

The sample problems have been selected to illustrate 
various features of NEC-BSC. Each problem has been 

selected in order to show specific code capabilities and 
convenient options, as well as some unresolved problems 
still present in the most recent version of the code. The 
collection of results provides an insight into the 
tradeoffs of the NEC-BSC use for indoor propagation 
prediction. 

A. COMPARISON OF NEC-BSC AND "PATCH" CODE RCS PREDICTIONS 

The problem of EM scattering by a smooth, perfectly 

conducting convex surface, is investigated in this section. 

The scattered field results obtained from the "PATCH" code 

[9] will be compared to the scattered field calculated by 

the NEC-BSC for the same geometry. The first example 

involves an infinitely thin electric dipole, shown as a red 

dot in Figure 6, located parallel and relatively close to a 

metallic circular cylinder. The cylinder axis coincides 

with the z-axis of the coordinate system. The operating 

33 




frequency is 2.54 GHz, which corresponds to the operating 
wavelength of 11.81 cm. The obstacle is a metallic 
cylinder 1 meter in length and having a 17.5 cm radius. 
The antenna is 1 meter from the surface of the cylinder. 
The observation points of interest form a horizontal circle 
of 360 points, 5 meters away from the cylinder axis, as 
shown in the Figure 6. 



Figure 6. Geometry of the PEC Convex Surface and the 
Dipole. 

The validity of NEC-BSC results is shown by comparing 
the NEC-BSC results with the solution obtained for the same 
problem by using the "PATCH" code. The results from the 
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two codes are shown combined in a single graph in Figure 7. 


Electric field measured by patch (blue) and NEC-BSC 



Figure 7. Comparison of Electric Fields Calculated by NEC- 
BSC and PATCH Codes. 

It can be clearly stated that the match is not 
perfect. Although the number of minima and maxima are the 
same and the "shape" of the field variation with angle is 
almost identical, there is approximately 1.7 dB difference 
between the values at each point of the graph. Moreover, 
in the very deep shadow area exactly behind the cylinder, 
there seems to be some differences as well, although at a 
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relatively low level. Despite the small differences, there 
is good enough agreement between the two to validate the 
NEC-BSC code results. 

NEC-BSC has been also extensively validated by its 
author. User's manual [34] contains several examples 
involving metallic cylinders, where patterns were taken in 
many different cuts, for many source positions and for 
different polarization. For these examples, comparison of 
NEC-BSC predictions with the measured results shows a very 
good agreement. 

B. SIMPLE ROOM WITH METALLIC WALLS 

This example considers a simple room with PEC walls 
and the following dimensions: 

• Length of six meters 

• Width of twelve meters 

• Height of twelve meters. 

The room walls are modeled using six PEC plates. 
Inside, there is a half-wavelength antenna located at 
(3,1,11) with respect to the origin, as shown in Figure 8 
(the antenna is depicted as a red dot) . The operating 
frequency is 2.4 GHz, which corresponds to a wavelength of 
12.5 cm. A large metallic cylinder, representative of a 
missile tube, is positioned in the room such that the 
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antenna illuminates only about half of the surface of the 
cylinder. Our objective, in this specific example, is to 
visualize the propagation paths, especially the non-LOS 
ones. The observation points are thus positioned in a 
semicircle and appear as gray dots in the figure below. 



Figure 8. Room With PEC Walls and a Single PEC Cylinder. 

Running the NEC-BSC program for even this relatively 
simple geometry results in a plethora of rays. Hence, even 
if we rotate or zoom into the figure, we will not get a 
clear picture of the individual propagation mechanisms 
involved. 
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NEC-BSC, however, provides the user with two very 
convenient options. The first is that the user can adjust 
the dynamic range of the results (in the ray menu) such 
that only the strongest rays are displayed (e.g., for LOS 
paths and single reflections only). For instance, in 
Figure 9 only the rays within 1 dB of the strongest ray are 
shown. Second, we can select a specific ray and transfer 
it to a figure, which only shows the geometry of the room 
(the "drag and drop" method). Therefore a "hand-picked" 
number of rays can be superimposed to the room geometry. 
To illustrate this for the room with the single cylinder 
example, only the diffracted rays to some certain 
observation points are shown in Figure 10. 



Figure 9. Propagation Paths for Signals Within 1 dB of the 
Strongest Signal Path. 
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Figure 10. Sample Diffracted Rays, "Isolated" from Figure 
7. 

C. COMPLEX ROOM WITH METALLIC WALLS 

Let us now consider a more complex room with PEC walls 
and the same dimensions as the room from the preceding 
section: 

• Length of six meters 

• Width of twelve meters 

• Height of twelve meters. 

Six PEC plates form the room. Inside, there is a 
half-wavelength antenna located at {3,3,1} with respect to 
the origin, depicted as the red dot. The operating 
frequency is 2.45 GHz, which corresponds to the wavelength 
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of 12.245 cm. Two large, vertical metallic cylinders, 
representing missile tubes, are placed in the center of the 
room. In addition, two large metallic tanks are located 
along the front and the rear wall of the room, as shown in 
Figure 11. The objective, in this specific example, is to 
obtain relative signal strength "map" of the room, 
revealing the areas where the received signal is strong, 
weak or virtually non-existent. 



Figure 11. Room With PEC Wall, Two PEC Cylinders and Two 
PEC Storage Tanks 







NEC-3SC does not have the 3D "relief map" plot 
capability to represent the calculated fields. A program 
has been written in Mathcad to create the relative signal 
strength "relief" (3D) map. More specifically, using the 
Notepad application, we convert the data file that NEC-BSC 
creates for the electric field at the observation points to 
an ASCII text file (*. txt extension). The ASCII data file 
is next imported into a Mathcad script, and the electric 
field magnitude is plotted in dB, as shown in the next 
Figure. 

The observation points are all in the same plane so 
that the two horizontal axes are the coordinates in the 
observation plane and the color-coded vertical axis (the 
"height") indicates the electric field strength. 



Figure 12. Signal Coverage of the Room, in Which Red 
Corresponds to 0 dB and Blue to -25 d3. 




A 


Figure 13. Top View of the Signal Coverage of the Room, in 
Which Red Corresponds to 0 dB and Blue to -25 dB. 

The red color identifies the areas of the highest 
signal strength while blue identifies the areas of weakest 
signal. It is important to remember that although the NEC- 
BSC will report signal strengths with very high dynamic 
ranges, the NEC-BSC user's manual [34] cautions that only 
about 30 dB of the dynamic range (30 d3 down from the 
highest signal level) should be considered reliable. 
Therefore, one can opt to place a threshold at 30, dB below 
the largest signal strength and thus establish a scale with 
red color indicating the largest signal strength and blue 
color indicating signal strengths at 30 dB lower (than the 
maximum). 
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As expected, the signal strength for observation 
points within PEC. objects (the two cylinders and the tanks) 
are represented by blue color, while the area close to the 
antenna has a large number of red "spikes" corresponding to 
the standing wave peaks. The pronounced standing waves are 
a consequence of strong reflection off the PEC walls. 


cyii 

.nders and tanks. 

The field 

strength maps 

also reveal 

the 

presence of the 

signal in 

the "shadow" 

zones area 

behi 

nd the cylinders, 

due to 

diffracted and 

"creeping" 


waves. 

D. OPENINGS 

In this section the performance of NEC-BSC is 
evaluated for openings such as doors, hatches, and windows. 
The opening will be constructed according to the method 
recommended in the NEC-BSC user's manual. An electrically 
large opening, showing in Figure 14, will be constructed 
using a number PEC plates attached together (by the virtue 
of shared corner pairs). These individual places form a 
larger flat plate with the opening in the middle. A half- 
wavelength antenna is positioned on one side of the opening 
and a set of observation points is defined on the opposite 
side as shown in Figure 15. In addition to the direct rays 
passing through the opening we expect to see the rays 
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Observation 



Figure 15- The Flat PEC Plate With the Opening in the 
Middle, as Constructed by NEC-BSC. 

The NEC-BSC results for the PEC plate with the opening 
show that, in addition to the rays diffracted from the 
edges of the opening, there are non-physical rays 
diffracted by the seams of the individual plates, as shown 
in Figures 16 and 17. 
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Figure 16. Sample Propagation Paths for the PEC Plate With 
the Opening. 



Figure 17. Non-Physical Rays "Diffracted" by the Plate 
"Seams 

Unfortunately, this problem was not resolved either 

by constructing the large plate using overlapping (by 0.1 

of the wavelength) individual plates, or by using the 

multiple olate {PM) command. We have therefore concluded 
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that the most likely cause is an unreported (and perhaps an 
undetected) "bug" in zhe code that: could afreet the 
accuracy cf the code in certain cases (for example for 


room-to-room 

propagation 

simulation) . A 'way to 

circumvent 

the "leakage 

through the 

seams" problem has been 

discovered 


by trial and error. 

In this approach one starts with a single, solid PEC 
plate. Next, a second plate is defined such that it 
coincides with the desired opening in the PEC plate. This 
second plate is assigned electrical properties 
(permittivity and permeability) of free space. The NEC-BSC 
apparently treats the second plate as an "opening," with 
direct and diffracted rays, but without the non-physical 
rays diffracted from the "seams" (because there are no 
seams) . Figure 18 shows the LOS and the diffracted paths 
for the "seamless" PEC plate with the opening. The 
difractions from the edges of the opening are evident. 
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igure 18. Propagation Paths for the "Seamless" Plate With 
he "Opening" Formed by a "Free Space Plate". 
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V. 


TYPICAL MISSILE ROOM 


A. GENERAL DERCRIPTION 

As mentioned, providing good signal coverage for 
wireless communications on a military ship can be a very 
challenging task. One of the objectives is to provide the 
best coverage with the smallest number of "access point" 
stations ("base" stations). To accomplish this, computer 
simulations may be used for the particular shipboard 
compartments to determine the signal coverage for various 
combinations of the numbers and locations of antennas (with 
multiple antennas it is assumed that they transmit 
identical signals; that is, the "multicast" operation is 
assumed). The requirement to minimize the computational 
time must be balanced against the desired accuracy of the 
solution. Therefore, the optimum level of geometric detail 
in the compartment models results from a tradeoff between 
the desired accuracy and the computational resources 
available (CPU time and memory). 

Simulation of a "typical" shipboard missile room is 
presented in this chapter. Although the geometry of missile 
room does not represent the missile compartment of any 
particular ship, it is typical of missile compartments for 
submarines. 
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The following software versions were used for the 
simulations: 

• NEC-BSC graphical user interface (GUI) called 
Workbench, version 4.1.37 and 

• NEC-BSC executable version 4.2-06. 

These were the most recent releases by Prof. Ron Marhefka, 
Ohio State University, at the time of completion of this 
thesis. 

In shipboard environments the majority of materials 
involved are metallic. Therefore, although lossy 

dielectric materials - are also present in shipboard 
compartments (people for example), propagation within 
shipboard compartments will be dominated by reflections and 
diffraction of metallic surfaces. Furthermore, the metals 
involved may be approximated as PEC. 

The large, parallelogram-shaped compartment shown in 
Figure 9 has the following dimensions: 

• Length of twenty meters 

• Width of ten meters 

• Height of ten meters. 

Six PEC plates form the compartment. Located at the first 
floor is a half-wavelength antenna, depicted as the red dot 
in the figures. The operating frequency is 2.4 GHz, which 
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corresponds to the wavelength of 12.5 cm. Six large, 
vertical metallic cylinders, representing missile tubes, 
are placed in the room. Each tube is 2 meters in diameter 
and 10 meters high. 

Furthermore, there are three large tanks, which cover 
the walls parallel to the x-axis. The front and the rear 
tanks' height are 2.5 and 0.5 meters respectively, while 
the middle one's width is 1 meter and the height is 5 
meters (same as the height of the first floor's ceiling). 
In addition to these six tanks, there is one large tank of 

2.5 meters height along the wall parallel to y-axis on the 
rear side of the room. 

One hatch is located in the middle of the front 
bulkhead. Since, hatches do not start at the floor level 
for the safety reasons (water containment), the hatch is 
positioned 30 centimeters from the floor and its height is 

1.5 meters. Two similar hatches are located on the first 
and the second overheads. 

The hatches can be modeled as either closed or opened. 
Therefore, propagation from one compartment section to the 
others through open hatches can be simulated. Finally, a 
nine-siep stairway, complete with metallic handrails, 
connects the first and the second floor. The geometry of 
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the NEC-BSC missile compartment model is shown in Figure 
19. The top view of the geometry of the compartment is 
provided in Figure 20. 



Figure 19. General View of a Typical Missile Shipboard 
Compartment. 


54 
































B. NEC-BSC MODEL SPECIFICS 

The NEC-BSC model described in the preceding section 
met all the requirements with respect to the limitations of 
the code, stated in Chapter III. In addition, all plates 
were modeled in a way that the code could infer which side 
of the plate was illuminated by the antenna. This was 
accomplished by defining all plates according to the right- 
hand rule described in Chapter V of the NEC-BSC manual. To 
eliminate NEC-BSC run-time "warnings" concerning the wall- 
attached plates of the rear tank, we had to define them as 
having corners not exactly identical, but positioned a 
small amount (0.5A) through the surfaces to which ' there 
were being attached. 

In the cases where we had to place identical objects 
in different locations, the NEC-BSC "Relative Rotate- 
Translate" (RR) command was used. This command enables the 
user to translate and/or rotate the coordinate system used 
to define the input data in order to simplify the 
specification of the set of missile tubes and the numerous 
stairs. Unlike the RT command, which always referred back 
to the reference coordinate system, the RR command 
accumulated the translations and/or rotations necessary for 
the desired positions. Therefore, at the end of the 
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procedure it was necessary to refer to the global 
coordinates again, by placing and then commenting out a 
test object at the origin that verified the accuracy of the 
procedure. 

Finally, this NEC-BSC model did not account for the 
presence of the people. Although presence and/or motion of 
people can cause up to 10 dB attenuation of the average 
signal strength [35], the trend towards reduced ship 
manning [36] through increased automation and sensor 
monitoring via ship-wide networks, means that the presence 
of a large number of personnel in a missile compartment is 
unlikely. 

C. SIMULATION RESULTS 

1. Antenna Located in the Middle of the Room 

In this case, both antenna and observation area are 
located on the first floor. The source is located in the 
middle of the compartment, 2 meters above the floor level, 
while the field observation points form a horizontal 
rectangle of 400 by 200 points (for a total of 80,000 
observation points) at the same height as the antenna. 
Recall that the operating wavelength is 12.5 cm, while the 
x and y dimensions of the compartment are 20 by 10 meters 
respectively. Therefore, the minimum spatial sampling that 
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allows us to capture wave minima and maxima requires a 
sampling step of one half of the operating wavelength, 
which corresponds to at least 320*160 = 51,200 observation 
points. Increasing the density of the observation points 
results in a more accurate "map" of the field, but at the 
expense of increased computational time. It was decided 
that 80,000 points represent a good compromise between the 
detail level and the computational time required for the 
solution. 
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Antenna’s 

location 



Figure 20. Top View of the Compartment. 

Despite the reasonable computing resources used (a PC 
with a single Pentium III 600 MHz processor and 768 
Megabytes of RAM), the simulation run time was still 
approximately seven hours for the 80,000 observation 
points. 
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By extracting the numerical results from. NEC-BSC and 
importing them to Mathcad, a "map" of the signal strength 
within the compartment was obtained, shown in Figure 21. 
It should also be noted that the threshold for the signal 
strength at the observation points was set at 40 dB below 
the peak observed signal strength. As stated before, the 
dynamic range for which the NEC-BSC has good accuracy is 
approximately 25 dB. Nevertheless, a larger dynamic range 
has been selected, in order for the shadow areas and the 
obstacles to be clearly visible in the signal strength 
"map." 



Figure 21. RF Energy Propagation Inside the Room {2 Meters 
Above the Floor). 
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2. Results for the Antenna in an Arbitrary Location 

Again, both antenna and observation points are located 
on the first floor. The antenna, however, is now placed at 
an "arbitrary" location (no longer in the middle of the 
compartment) 1 meter above the floor, while the observation 
points form a horizontal rectangle of 400 by 200 points 
(for a total of 80,000 observation points) and 1.5 meters 
above the floor. The top view of the compartment is 
provided in Figure 22. The arrow indicates the antenna 
location. 
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Antenna’s 

location 



1-1 

Figure 22. Top View of the Missile Compartment With the 
Antenna in an Arbitrary Location. 

The "maps" of the signal strength are shown in Figures 
23 and 24, for the dynamic ranges of 33 and 22 dB 
respectively. In the 33-dB dynamic range "map" we notice 
the apparent signal "leakage" into the missile tubes and 
the tanks. Since the tubes and the tanks are non- 
penetrable (modeled as PEC) these are clearly consequence 
of the code's limited dynamic range. Reducing the dynamic 


61 




range to 22 dB eliminates most of the "non-physical" rays, 
as evident in Figure 24. 



Figure 23. Signal Strength "Map" for the Missile Room Using 
a Dynamic Range of 33 dB. 



Figure 24. Signal Strength "Map" for the Missile Room Using 
a Dynamic Range of 22 dB. 
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3. 


Results for two Antennas in the Missile Room 


Next, the results are next presented for a "typical" 
wireless LAN installation where a number of antennas may be 
used to increase the uniformity of the signal coverage 
within the user space. In particular, it is assumed that 
the two antennas transmit identical signals, a situation 
commonly referred to as "multicasting." The two antennas 
are centrally located, close to the ceiling, at the height 
of 4 m above the floor, and 8 m apart from each other, as 
shown in Figure 25. The 80,000 observation points again 
form a horizontal rectangle (400 by 200 points) at 1 meter 
above the floor (approximately waist level). 
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Position of the two antennas 



I-1 

Figure 25. Top View of the Compartment With the Two 
Antennas. 

By extracting the NEC-BSC solution (numerical results 
for the field strength) and importing it into Mathcad, we 
obtained the "map" of the signal strength within the 
compartment, shown in Figure 26. The dynamic range for 
this "map" was limited to 25 dB. As expected, the "map" 
shows that a much more uniform coverage is obtained with 
two antennas. There is even coverage behind the 

cylinders, due to reflected and "creeping" wave 
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contributions. 


The computer running time for this 


simulation was approximately 12 hours. 



Figure 26. Signal Propagation Inside the Missile Room, 4 
Meters Above the Floor and Using a Dynamic Range of 25 dB. 

4. Floor-to-Floor Propagation 

The simulations have confirmed the 20 to 30 dB dynamic 
range limitation of for NEC-BSC stated in the NEC-BSC 
User's Manual. The limited dynamic range suggests that 
NEC-BSC will not be well suited for applications where the 
desired signal levels could be very low, such as in the 
case of floor-to-floor propagation. Nevertheless, 

simulations were performed to test the NEC-BSC performance 
for room-to-room and floor-to-floor propagation. Presented 
here are the NEC-BSC results for floor-to-floor propagation 
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signal strength for the floor-to-flcor propagation, shown 
in Figure 27. Since floor-to-floor propagation incurs a 
considerable loss (no LOS path) we had to use 75-dB dynamic 
range for the "map." This is clearly beyond the stated 
reliable dynamic range of the NEC-3SC and is manifested by 
what appears to be the signal "leakage" at the "seams" 
between the plates forming the fioor/ceiling and the 
"seams" between the tubes and the floor/ceiling, evident in 
Figure 28. Nevertheless, the code properly predicts the 
highest signal strength for the vicinity of the stairwell 
opening. The computer run time for this simulation was 
approximately 12 hours. 



Figure 28. Signal Propagation Inside the Second Floor, 
Caused by an Antenna Located in the First one and Using a 
Dynamic Range of 75 dB. 
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However, this simulation does give us some insight 
about the diffraction and reflection mechanisms that take 
place in communication between floors. These phenomena, 
which dominate in room-to-room cases, can be generally seen 
in Figure 28. The run time for this program was 
approximately 12 hours. 

D. PROPAGATION STATISTICS FOR THE MISSILE COMPARTMENT 

1. Theoretical Background 

Thus far, we have presented the signal strength "maps" 
for particular compartments using the numerical output 
results of the NEC-BSC and Mathcad. In this section the 
objective is to extract the statistics of the excess path 
gain/loss (relative to free space path loss) for the entire 
floor of a particular compartment and for different antenna 
placements. The result will be a random variable with a 
certain probability density function (pdf) expressing the 
likelihood of a particular excess path gain/loss. 

The path loss for free space can be easily calculated 
for any distance between the transmitter and an observation 
point. Adding a random variable with a known pdf to the 
deterministic free space loss allows us to calculate the 
mean value, dispersion, etc. for the compartment path loss 
random variable. This applies, in statistical sense, to 


68 



any two points (the transmitter and the observation point) 
inside the particular compartment, regardless of their 
exact location as long as the distance between them is the 
same. 

For each transmit antenna location, two NEC-BSC input 
files were created: one for the compartment and the other 
for the free space of exactly the same shape and volume as 
the compartment. All other parameters (antenna type and 
position, radiated frequency, observation points placement) 
were kept identical. The two output ASCII files, resulting 
from the execution of the input files, were formatted using 
the Windows WordPad editor and imported into Mathcad'. Two 
column matrices were created in Mathcad for each input 
file, with the first column containing the relative 
distance between the transmit antenna and the observation 
point and the second column containing the corresponding E- 
field strength. 

The matrices are then rearranged in the following 
manner. First, the distances (from minimum to maximum) are 
discretized into distance "bins," each bin 50 cm (4 
wavelengths) wide. Next, all values within the same 
distance bins are placed on the same row of a new "matrix" 
(this "matrix" has a variable number of columns). The 
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distance "bins" geometrically correspond to concentric 
spherical shells of 50 cm thickness, all centered at the 
transmit antenna location. The number of observation 
points within each "shell" varies with the shell distance 
from the "origin" (the transmit antenna location). 

The next step is to form the ratios of the compartment 
and "free space" fields at the same observation points. 
Expressing this ratio in dB creates a new variable that is 
referred to as the "excess compartment gain" R. Positive 
values of the excess compartment gain indicate that the 
signal in the compartment exceeds the free space signal at 
the same location, while the negative values of the excess 
compartment gain indicate that the signal in the 
compartment is weaker than the free space signal at the 
same location. 


jR = 20 log 


'room 

Efree 


(5.1) 


In Equation 5.1, R denotes the "excess room gain" with 
respect to free space. Since for a given distance bin 
there is a large number of values (because of a typically 
very large number of observation points) the values of R 
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from any given bin form a random variable. The random 
variables for the individual distance bins can be further 
combined to form a single random variable that now 
incorporates the excess gain statistics for the entire 
compartment, regardless of the distance. 

The probability density function (pdf) and the 
cumulative distribution function (cdf) for the excess gain 
random variable are compared to the Gaussian pdf and cdf 
with the same mean and variance. The pdf's and the cdf's 
are plotted versus the excess gain in dB. The 

distributions obtained via simulation are plotted on the 
same graphs with the Gaussian distributions, to illustrate 
the "match" between the two. In addition, the histograms 
for the distance groups are also presented. 

In order to obtain the aforementioned plots, we first 
had to resolve an important issue. Invoking the volumetric 
plot command (VN) results in a number of observation 
points. The observation points are uniformly spaced and 
thus would be positioned inside as well as outside non- 
penetrable objects like the cylinders, the tanks, even the 
steps of the stairways. Since the interiors of non- 
penetrable objects are not considered as "legitimate" 
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transceiver locations, the observation points within non- 
penetrable objects had to be eliminated. 

A simple attempt to circumvent this problem by the 
multiple selective use of the same command failed since the 
code executed only the last of the commands. Therefore, we 
had to apply a rather time consuming procedure: we 
truncated the missile room and the free space area into a 
number of obstacle-free observation regions. We ran the 
program once for each of the observation regions, and then 
"stacked" (using Mathcad) the data resulting from each run 
to form a single "output" file. 

The results of the simulations presented in this 
chapter indicate that the Gaussian distribution best 
describes the missile compartment excess gain (in dB) 
random variable. It has been reported [1] that the 
lognormal distribution also applies to indoor propagation 
between observation points that have different distances 
from the RF transmitter. 

Furthermore, it has been reported recently [38] that 
the lognormal distribution applies to the "same distance 
points" that form concentric cycles around an arbitrarily 
placed source for indoor propagation within a building. 
The same approach has been applied in this, thesis to 



metallic compartments/ with similar results, that is the 
lognormal distribution for the excess gain random variable. 
Recall that the missile compartment is composed of PEC 
materials, hence it can be characterized as highly 
reflecting/diffracting environment. 

Due to the absence of absorption and refraction, 
reflected, diffracted, and "creeping" waves dominate in the 
propagation mechanisms. Each of the multiple reflections 
that a single "ray" experiences within the compartment 
contributes a factor (the reflection coefficient to the 
magnitude of the ray at the particular observation point). 
Therefore the observed ray magnitude is a product of a 
number of factors. The logarithm of the product (that is 
the path gain/loss in dB) is equal to the sum of the 
logarithms of each term in the product. Since each 
reflection is by and large independent from the others 
(involving different locations, objects, and incident 
angles) the sum of logarithms tends to be a normally 

distributed random variable, hence the lognormal 

distributions. 

2. Results 

Let us now examine the results for several 
transmitting antenna locations. The observation points are 


73 



always located on a horizontal plane, but the observation 
plane height {above the floor level) may differ for each 
individual case. 

2.1 First Location 

The transmitting antenna is located at (14,5,1) and 
the observation points lay in a horizontal plane at Z = 
1.5m, as show in Figure 29. The various cells are shown in 
Figures 30 through 40. 


Location of 
Antenna 



Figure 29. Side View of the Compartment, the Observation 
Area, Forming a Line, and the antenna at location 1. 
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Figure 30. Top View of the Compartment, the First 
Observation cell, and the antenna at location 1. 
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Figure 31. Top View of the Compartment, the Second 
Observation cell, and the antenna at location 1. 
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Figure 32. Top View of the Compartment, the Third 
Observation cell, and the antenna at location 1. 
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Figure 33. Top View of the Compartment, the Fourth 
Observation cell, and the antenna at location 1. 
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Figure 34. Top View of the Compartment, the Fifth 
Observation cell, and the antenna at location 1. 
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Figure 35. Top View of the Compartment, the Sixth 
Observation cell, and the antenna at location 1. 
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Figure 36. Top View of the Compartment, the Seventh 
Observation cell, and the antenna at location 1. 
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Figure 37. Top View of the Compartment, the Eighth 
Observation cell, and the antenna at location 1. 
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Figure 38. Top View of the Compartment, the Ninth 
Observation cell, and the antenna at location 1. 
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Figure 39. Top View of the Compartment, the Tenth 
Observation cell, and the antenna at location 1. 
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Figure 40. Top View of the Compartment, the Eleventh 
Observation cell, and the antenna at location 1. 

The excess gain (in dB) pdf and cdf, obtained by 
averaging over all distance groups, are shown in Figures 41 
and 42, respectively. Plotted on the same graphs are the 
lognormal (Gaussian) pdf and cdf, to visualize the 
agreement between the lognormal distribution model and the 
data obtained by simulation. 
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Antenna at location 1 



Room gain (dB) 

_average pdf 

. .. gaussian 


Figure 41. Pdf of the Excess Gain and the Lognormal 
(Gaussian) Pdf, Transmitting Antenna at Location 1 








Antenna at location 



Figure 42. Cdf of the Excess Gain and Lognormal (Gaussian) 
Cdf, Antenna at Location 1 

As shown, there is a good match between the excess 
gain pdf (obtained by simulation) and the lognormal pdf. 
The oscillations ("dispersion") of the simulation results 
about lognormal distribution are inherent to the 
statistical nature of the simulation results. The 
"dispersion" can be reduced by increasing the number of 
observation points but at the expense of increasing the 
computational times. The agreement is even better for the 
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cdf's, which could have been expected as the integration 
"averages out" the statistical dispersion of the pdf- Note 
that the pdf oscillations are, on the average, equally 
likely to be either above or below the lognormal pdf, hence 
they tend to cancel out when the pdf is integrated to 
obtain the cdf. From the cdf we note that the compartment 
(for this particular antenna location) has a 0.5 
probability of behaving exactly like free space; that is, 
it experiences neither gain nor loss relative to free space 
propagation. 

2.2 Second Location 

In this case the transmitting antenna is located at 
(2,2,2) and the observation points are in a horizontal 
plane at Z = 2m. The results are illustrated in the Figures 
4 3 and 4 4. 
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this "thresholding" causes a distinct peak to appear in the 
pdf, distorting the pdf's shape relative to the expected 
lognormal distribution. To illustrate the effects of the 
dynamic range limiting on the pdf obtained via simulation, 
the pdf for the 40-dB dynamic range is shown in Figure 45. 
Compared with the 30-dB dynamic range pdf, the 40-dB 
dynamic range pdf shows a more prominent peak, shifted 
further to the left (larger values of excess path loss). 

In order to mitigate the effects of the dynamic range 
limitation inherent to the NEC-BSC (as a mater of fact all 
simulation codes have a finite dynamic range, due to the 
"numerical noise") a "fix" has been implemented in Mathcad. 
The essence of the "fix" was to replace the value of the 
"offending" peak with the value that corresponds to the 
lognormal distribution and then scale the resulting 
distribution such that its integral is equal to 1 (as it 
has to be for any pdf). Shown in Figure 45 are the pdf's 
for the 30 and 40 dB dynamic ranges and the "corrected" (in 
the manner described above) pdf for the 30-dB dynamic 
range. It is evident that the outlined procedure results 
in a pdf that is again close to the lognormal pdf. 
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2.3 Third Location 


In this case the transmitting antenna is located at 
(2, 2 , 2.2) and the observation points are in a horizontal 
plane at Z = 2.2m (at the same height as the transmitting 
antenna) . Note that the difference between this and the 
preceding case is that the antenna and the observation 
points have been translated upward by 20 cm. The reason to 
do this is to examine the consequences of a translation by a 
relatively small distance, commensurate with the operating 
wavelength (approximately 1.5 wavelengths in this case). In 
addition, the observation points do not fall on the wide 
horizontal steps of the stairway. The results are shown in 
Figures 46 and 47. 
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Figure 46. Pdf of the Excess Gain and the Lognormal 
(Gaussian) Pdf, Transmitting Antenna at Location 3. 
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Figure 47. Cdf of the Excess Gain and the Lognormal 
(Gaussian) Cdf, Antenna at Location 3. 






















2.4 Fourth Location 


In this case the transmitting antenna is located at 
(4,3,1) and the observation points are in a horizontal plane 
at Z = 1 m (at the same height as the transmitting antenna). 
This is the last case considered, for which the transmitting 
antenna is placed at an "arbitrary" location (not in a plane 
of symmetry, or the like). However, the observation points 
are still in a horizontal plane, at the same height as the 
transmitting antenna. Following the same procedure as for 
the preceding three locations, we obtain the results shown 
in Figures 50 and 51. 
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Figure 50. Pdf of the Excess Gain and the Lognormal 
(Gaussian) Pdf, Transmitting Antenna at Location 4. 





















Average cdf 
Gaussian cdf 

• • * 0.4 




Room gain (dB) 


25 30 


Figure 51- Cdf of the Excess Gain and the Lognormal 
(Gaussian) Cdf, Antenna at Location 4. 


The data for the four antenna locations used in the 
simulations are presented in Table 5.1. 


Antenna 

Antenna 

| Observation 

j 

location 

Coordinates ; 

! points 


Excess 


Excess gam 



Second 


Third 


2 , 2,2 


1.5 


Z=2 


gain St. Deviation 



Mean (dB) 


- 0.10 


1.43 


2.91 


0.5 



(dB) 


5.31 


4.85 


3.95 


3.9 


’able 5.1. Data for the Four Propagation Simulations. 
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VI. UNRESOLVED ISSUES, RECOMMENDATIONS, AND 

CONCLUSIONS 

A. UNRESOLVED PROBLEMS 

1. Floor-to-Floor Propagation 

The case in which the antenna was located inside the 
first floor of the missile room, and the observation points 
in the second (and vice versa), could not be thoroughly 
examined due to the problems caused by the "ray leakage," 
as reported in Chapter IV. Therefore the compartment-to- 
compartment propagation statistics were not addressed. 
Note that the limited dynamic range is also an issue with 
regard to compartment-to-compartment coupling that "by 
definition" involves large attenuation. 

2. Problematic PN Command 

The Near Zone Pattern (PN) command enabled us to 
define the near zone pattern coordinate system, the 
volumetric pattern cut, and the spatial range of desired 
observation points. However, as reported in Chapter V, the 
multiple and selective use of the PN command to exclude the 
undesired observation points (inside non-penetrable 
objects) was unsuccessful, because the code executed only 
the last of the multiple PN commands. The truncation of 
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the room into several regions, exclusive, of non-penetrable 
objects, did work but at the expense of excluding some 
points located near the non-penetrable obstacles and it 
increased time requirements for processing of multiple data 
files. 

3. Running Time 

All simulation programs executed in Chapter V were 
time consuming and had run times, ranging from seven to 
twelve hours. As expected, time requirements depended on 
the complexity of the compartment, the number of the 
observation points and the number of antennas used. 

B. RECOMMENDATIONS 

All recommendations for further study are made in the 
light of existing information, based on NEC-BSC Workbench 
version 4.1.37 and executable NEC-BSC FORTRAN program 
versions 4.2-06. 

Some commands, namely BF, FT, LY, NC, NG, NM, NP, NR, 
NS, NX, PT, RI, SI, TA, TD, TW, TY, VM, were not supported 
by the Command Wizard window of the code's Graphical User 
Interface (GUI). Moreover, the function of the command ZC, 
provided by Workbench, remains mysterious, since there is no 
description for it, either "on-line" or in the User's 
manual. 
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In addition, the Zoom In command in the main toolbar 
was inconvenient, because the user was not able to specify 
the desired scale, which by default was inadequately large. 
The accuracy of the Zoom was also questionable, since it 
altered the location of certain objects, presented the 
cylinder caps as polygons, and showed some rays diffracted 
from cylinders as if the rays penetrated the cylinders and 
then diffracted from the inside of the cylinders. 

It should also be stated that the very useful and 
convenient "drag and drop" option couldn't be used in all 
cases; the user could not always distinguish and select an 
individual ray among a plethora of others. Finally, the 
coordinate system of the code, which by default placed the 
origin in the left upper corner of the display, used 
unknown units and was difficult to work with. 

A recent, extensive study [38] includes a number of 
recommendations that would significantly improve the 
applicability of the NSC-BSC for indoor propagation 
simulation. In addition to these, parallel execution 
option [37] would also decrease the run-times, especially 
in complex cases that approach the limitations of the code, 
as reported in Chapter III. 
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Moreover, a considerable improvement in the 
calculation speed and memory requirements could be made if 
the near-field computation by NEC-BSC did not include both 
the Electric and the Magnetic field components 
simultaneously. In situations where the computation of the 
signal strength was the goal, the calculation of the H- 
field could not be skipped. The user should be able to 
choose the particular quantity that the code should 
calculate in order to reduce the computational time. 

Another factor that would improve code performance is 
the support of all commands in the Workbench. In addition, 
the Users Manual should be completed and should include all 
the figures omitted in the previous versions, along with 
more examples that illustrate the workflow. 

Finally, the creation of the problem geometry will be 
greatly simplified if future versions of NEC-BSC would 
support data import in a standard CAD file format {.dxf, 
.dwg, etc.). 

C. CONCLUSIONS 

The primary objective of this research was to provide 
simple statistical models for indoor propagation in complex 
environments, such as shipboard compartments. The emphasis 
in this research was on the basic principles and a simple 
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statistical description for a representative shipboard 
environment rather than on encyclopedic compilation for 
indoor propagation data for vastly diverse complex 
environments. Previous work using NEC-BSC for indoor 
propagation simulation, involving environments with no 
material uniformity [38] showed satisfactory agreement with 
both the results reported in [40], and with the 
measurements conducted in the Bullard Hall at the Naval 
Postgraduate School in Monterey, California [38]. 

In our case, as well as for the Bullard Hall case 
[38], the "room excess gain" random variable was defined to 
describe how the propagation in the indoor environment 
differs from the deterministic (formula-predictable) free 
space propagation. The simulation results in both cases 
show that the excess gain (in dB) random variable could be 
approximated as lognormal (Gaussian) distributed [39]. 
Note that this random variable is representative for all 
distance combinations and locations of the source and 
observation points within a single compartment. Although 
the shape of the probability density function for the 
excess gain (in dB) random variable does not depend on the 
particular choice of the source location or the compartment 
details, the parameters defining the pdf, such as the mean 
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value and the variance, do depend on the position of the 
source and the observation points relative to the obstacles 
present in the compartment. 

It appears that the key parameter in defining the mean 
value of the excess gain is the relative number of 
observation points that are shadowed by the obstacles and 
thus experience non-line-of-sight propagation attenuation. 
If there is a large number of observation points in deep 
shadows, which depends primarily on the location of the 
source relative to the obstacles in the compartment, the 
mean value of the excess gain is negative, indicating a 
loss (on the average) relative to the free space 
propagation. If on the other case most of the observation 
points are either direct line-of-sight or are obstructed by 
low-loss penetrable obstacles, the mean value of the excess 
gain (in dB) becomes positive, indicating that the signal 
strength in that compartment exceeds, on average, signal 
strength in free space. 

The later can be physically understood by observing 
that the signal energy in free space propagates outward, 
with the corresponding decrease in the power density, while 
in an indoor environment - some of the energy is reflected 
back by the walls, window panes, etc., thus increasing the 
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average power density within the compartment- The 
particular compartment analyzed in this research and for 
the several transmitting antenna locations is by and large 
"free space" neutral, with the excess gain mean value 
between -0.1 and +2-9 dB, depending on the source location. 
The restrictions of time and computer resources limited 
this research to only one complex compartment. 
Furthermore, the computer resources available would not 
allow the placement of observation points at separations 
below "two per operating wavelength" to allow very fine 
sampling of the fields inside the compartments. 

These facts, as well as the code limitation with 
respect to the dynamic range, should be taken into account, 
and offer plenty of opportunities for further research. 
Nonetheless, the results definitely confirmed the 
applicability of the excess gain random variable approach in 
general and the lognormal distribution model in particular. 
NEC-BSC is generally accepted as a very useful tool for 
modeling a variety of radiation and scattering phenomena. 
Over the years it was reported to provide reliable results 
especially when used for outdoor environments. Recently, 
other studies ([38] and the present one) showed quite 
promising results in using NSC-BSC for modeling indoor 
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propagation as well. However, definite conclusions require 
extensive further investigation and modeling of a multitude 
of diverse indoor environments. 
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APPENDIX A. MATIiAB FILES 

This is the Matlab file used to read the data, resulted 
from "PATCH" software and create the plots necessary for 
the comparison between that and NEC-BSC (Case: cylinder and 
dipole) 

% plot near field data 
clear 

rad=pi/180; 
iff=l; 
if iff==0 
load amptffl.m 
load amppffl.m 

arrrptff=amptf fl; amppff=amppffl; 

P=[0:180] ; 
figure(1),elf 

plot(P,amptff(1:181)/P,amppff(1:181),' —'),grid 
xlabel(’Phi, deg') 

ylabel('Total electric field (dB)’) 

title('Far field, a=0.175m, L=lm, theta=90, dip @ phi=0') 

axis{[0,180,-30,0]) 

end 


% currents (get from outnear files) 

Js=0.17859E+01-rj*0.20941E-f01; iplt=0; % 2.54 GHz outnear 

2 

% Js=0.0041847-fj *0.40483E+00; iplt=0; % 1 GHz outnear 3 
% Js=0.18605E-f01-j*0.83378; iplt=l; ' % 1 GHz dipole 
outnear 6 (20 segs) 

% Js=0.96774E+01-j*0.48566E+01; iplt=l; % 1 GHz skinny 

dipole 20 segs outnear 7 

% Js=2.4268-j*0.62429; iplt=0; % outnear 9 (longer dip 
for 1 GHz) 

I=Js*.005; % I=Js*(width of dipole) 

V=l; 

Z=V/I; 

% D=20*logl0(abs(I)*.19); % for outnear 3 L is not wavl/2 

but 0.19 

D=20*logl0(abs(I));% factor of nonunit amplitudecurrent 
load nearpts2.m 
nearpts=nearpts2; 

[nrow,ncol]=size(nearpts); 
npts=nrow/2; 
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Ez (1:npts)=nearpts(1:2:nrow,7)+j*nearpts(1:2:nrow, 8); 

Ey{1:npts)=nearpts(1:2:nrow,5)+j*nearpts(1:2:nrow,6); 

Ex(1:npts)=nearpts(1:2:nrow,3)+j*nearpts(1:2:nrow,4); 

if iplt==0, 

load nearfldptsl80.m; nearfldpts=nearfldptsl80; 
end 

if iplt==l, load nearfIdptsz.m; nearfldpts=nearfldptsz; 
end 

x=nearfIdpts(:,1); y=nearfldpts(:,2); z=nearfIdpts(:,3); % 

set for 180 points 
R=sqrt(x. A 2+y. A 2+z. A 2); 

T=acos(z./R); 

P=atan2 (y,x) ; 

U=cos(T).*cos(P); 

V=cos(T).*sin(P); 

W=-sin(T); 

Y=cos(P); 

Z=sin(P); 

Et=U. * transpose (Ex) -fV. * transpose (Ey) +W. * transpose (Ez) ; 

Ep=-Z.*transpose(Ex)+Y.*transpose(Ey); 

Emax=sqrt(Et.*conj(Et)); % +Ep.*conj(Ep)); 

% Smax=sqrt (Ex. *conj (Ex) -rEy. *conj (Ey) 4-Ez . *conj (Ez)) ; 
Edb=20*logl0(Emax)-D; 
load gtddat.m 

E=gtddat(:,6); E(361)*E(1); 

if iplt==0 % plot cylinder and dipole 

E2=[Edb fliplr(Edb)]; % mirror image of pattern 

P2=[P/rad 360-fliplr(P/rad)]; 
figure(2),elf 

plot(P2,E2, 'b ', [0:360],E, 'r —'),grid 
xlabel('Phi, deg') 
ylabel('Total electric field (dB)') 

title('Electric field by patch (blue) and NEC-BSC (red)') 
end 


if iplt==l % just plot dipole data 

for i=l:length(T) 

TT(i)=T(i); 

if P(i)==pi, TT(i)=2*pi-T(i); 

end 

end 
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plot(TT/rad,Edb,'- 7 ) 7 hold on 
xlabel('Theta, deg') 
ylabel{'|Ej, dB') 

title{'half-wave dipole field at R=5m (f=l GHz)') 
axis([0,360,-20,30]) 
figure(3),elf 

polardb((TT/rad).',Edb,'-') 

title{'short dipole field at R=5m (f=l GHz)') 

% analytical dipole field (current is I from above) 
Et=377/2/pi/5*abs(cos(pi/2.*cos(TT))./sin(TT+le-3)) 
Etdb=20*logl0(Et); 
figure(2) 

plot (TT/rad, Etdb, ' —') 
axis([0,360,-20,30]) 
hold off 
end 

End of all Matlab files 
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APPENDIX B. NEC-BSC INPUT FILES 


CE: This is the input file used in the comparison between 
CE: NEC-BSC and PATCH software in the case of the cylinder 
CE: and the dipole 

ZD: 

0 

UN: 

1 

FR: 2.54 
2.54000 


US: wavelengths 
0 

SG: 


1 . 00000 , 

0 . 000000 , 

0.000000 


0 . 000000 , 

0 . 000000 , 

90.0000, 

0.000000 

-2, 0.500000, 0.000000 
1 . 00000 , 0.000000 


CG: 

0 . 000000 , 

0 . 000000 , 

0.000000 


0 . 000000 , 

0.175000, 

0 . 000000 , 

0.175000 

90.0000, 

0.000000 

-0.500000 

, 90.0000, 

0.500000, 

90.0000 

PN: 

0 . 000000 , 

0 . 000000 , 

0.000000 


0 . 000000 , 

F 

0 . 000000 , 

90.0000, 

0.000000 

5.00000, 

90.0000, 0 

.000000 


0 . 000000 , 

0 . 000000 , 

1.00000 



360 

PL: 

T 

LY 

T 

XQ: Execute the program 

EN: End 
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CE: This is the input file used in the simple metallic room 

CE: Six box room with a cylindrical obstacle inside it 

ZD: type of material: PEC 
0 

FR: we will use a frequency of 2.4 GK 
2.40000 

US: 

0 


SG: Source geometry (3,1,11). 

3.00000, 1.00000, 11.0000 
0.000000, 0.000000, 90.0000, 0.000000 
-2, 0.500000, 0.000000 
1.00000, 0.000000 


PG: First plate 
4, 0 . 

0 . 000000 , 0 . 000000 , 0.000000 
0 . 000000 , 12 . 0000 , 0.000000 
6 . 00000 , 12 . 0000 , 0.000000 
6 . 00000 , 0 . 000000 , 0.000000 

PG: second plate 
4, 0 

0 . 000000 , 0 . 000000 , 0.000000 
6.00000, 0.000000, 0.000000 
6 . 00000 , 0 . 000000 , 12.0000 
0 . 000000 , 0 . 000000 , 12.0000 

PG: Third plate 
4, 0 

6 . 00000 , 0 . 000000 , 12.0000 
6 . 00000 , 0 . 000000 , 0.000000 
6 . 00000 , 12 . 0000 , 0.000000 
6 . 00000 , 12 . 0000 , 12.0000 

PG: Forth plate ==front 
4 , 0 

6 . 00000 , 12 . 0000 , 12.0000 
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6 . 00000 , 12 . 0000 , 0.000000 
0 . 000000 , 12 . 0000 , 0.000000 
0 . 000000 , 12 . 0000 , 12.0000 

PG: Fifth plate 
4, 0 

6 . 00000 , 0 . 000000 , 12.0000 
6 . 00000 , 12 . 0000 , 12.0000 
0 . 000000 , 12 . 0000 , 12.0000 
0 . 000000 , 0 . 000000 , 12.0000 


PG: Sixth plate 
4, 0 

0 . 000000 , 0 . 000000 , 12.0000 
0 . 000000 , 12 . 0000 , 12.0000 
0 . 000000 , 12 . 0000 , 0.000000 
0 . 000000 , 0 . 000000 , 0.000000 

CG: Put a cylinder as an obstacle inside the room 
3.00000, 6.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 6.00000, 90.0000 


PN: Near zone cut 
3.00000, 6.00000, 3.00000 
0.000000, 0.000000, 90.0000, 0.000000 
? 

1.50000, 90.0000, 0.000000 
0 . 000000 , 0 . 000000 , 1.00000 
180 

PL: 

T 

LY 

T 

XQ: Execute the program 

EN: End 

CE: This is the input file used in the complex metallic 
room 

CM: Six box room (6,12,12) with two cylindrical obstacles 

CE: and two tanks inside 
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ZD: type of material: PEC 
0 

FR: we will use a frequency of 2.45 GH 
2.45000 

US: 

0 


SG: Source geometry (3,3,1). 

3.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-2, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT at Z=1.5 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

0.000000, 0.000000, 1.50000 
0 . 100000 , 0 . 100000 , 0.000000 
1, 60, 120 
PG: First plate 
4, 0 

0 . 000000 , 0 . 000000 , 0.000000 
0 . 000000 , 12 . 0000 , 0.000000 
6 . 00000 , 12 . 0000 , 0.000000 
6.00000, 0.000000, 0.000000 


PG: second plate 
4, 0 

0.000000, 0.000000, 0.000000 
6 . 00000 , 0 . 000000 , 0.000000 
6 . 00000 , 0 . 000000 , 12.0000 
0 . 000000 , 0 . 000000 , 12.0000 


PG: Third plate 
4, 0 

6 . 00000 , 0 . 000000 , 12.0000 
6 . 00000 , 0 . 000000 , 0.000000 
6 . 00000 , 12 . 0000 , 0.000000 
6 . 00000 , 12 . 0000 , 12.0000 


PG: Forth plate ==front 
4, 0 
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6 . 00000 , 12 . 0000 , 12.0000 

6 . 00000 , 12 . 0000 , 0.000000 

0 . 000000 , 12 . 0000 , 0.000000 

0 . 000000 , 12 . 0000 , 12.0000 


PG: Fifth plate 
4, 0 

6 . 00000 , 0 . 000000 , 12.0000 
6 . 00000 , 12 . 0000 , 12.0000 
0 . 000000 , 12 . 0000 , 12.0000 
0 . 000000 , 0 . 000000 , 12.0000 

PG: Sixth plate 
4, 0 

0 . 000000 , 0 . 000000 , 12.0000 
0 . 000000 , 12 . 0000 , 12.0000 
0 . 000000 , 12 . 0000 , 0.000000 
0 . 000000 , 0 . 000000 , 0.000000 

CG: Put a cylinder as an obstacle inside the room 
1.50000, 6.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
O'.SOOOOO, 0.500000 

0.000000, 90.0000, 12.0000, 90.0000 


CG: Put a second cylinder as an obstacle inside the room 
4.50000, 6.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
0.500000, 0.500000 

0.000000, 90.0000, 12.0000, 90.0000 


PG: first tank front 
4, 0 

2 . 00000 , 1 . 00000 , 0.000000 
2 . 00000 , 1 . 00000 , 6.00000 
4.00000, 1.00000, 6.00000 
4.00000, 1.00000, 0.000000 


PG: first tank side 
4, 0 

4.00000, 0.000000, 6.00000 
4 . 00000 , 0 . 000000 , 0.000000 
4.00000, 1.00000, 0.000000 
4.00000, 1.00000, 6.00000 

PG: first tank side 
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4 , 0 

2 . 00000 , 0 . 000000 , 6.00000 
2 . 00000 , 1 . 00000 , 6.00000 
2 . 00000 , 1 . 00000 , 0.000000 
2 . 00000 , 0 . 000000 , 0.000000 


PG: first tank up 
4 , 0 

2 . 00000 , 0 . 000000 , 6.00000 
4 . 00000 , 0 . 000000 , 6.00000 
4 . 00000 , 1 . 00000 , 6.00000 
2 . 00000 , 1 . 00000 , 6.00000 


PG: second tank front 
4 , 0 

2 . 00000 , 11 . 0000 , 0.000000 
4 . 00000 , 11 . 0000 , 0.000000 
4 . 00000 , 11 . 0000 , 6.00000 
2 . 00000 , 11 . 0000 , 6.00000 

PG: second tank side 
4 , 0 

4 . 00000 , 11 . 0000 , 6.00000 
4 . 00000 , 11 . 0000 , 0.000000 
4 . 00000 , 12 . 0000 , 0.000000 
4 . 00000 , 12 . 0000 , 6.00000 

PG: second tank side 
4 , 0 

2 . 00000 , 11 . 0000 , 6.00000 
2 . 00000 , 12 . 0000 , 6.00000 
2 . 00000 , 12 . 0000 , 0.000000 
2 . 00000 , 11 . 0000 , 0.000000 


PG: second tank ud 


4 , 0 

2 . 00000 , 

4 . 00000 , 

4 . 00000 , 

2 . 00000 , 

PL: 

T 

XQ: 


11 . 0000 , 

11 . 0000 , 

12 . 0000 , 

12 . 0000 , 


6.00000 

6.00000 

6.00000 

6.00000 
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EN: Ena 


CE: This is the input file used for openings 

CE: right hand ok, commands ok problem with penetration 

FR: We will use a frequency of 2.4 GHz (ISM band) 
2.40000 


US: 

0 


PG: platel up 
4 , 0 

6 . 00000 , - 6 . 00000 , 12.0000 
6 . 00000 , - 6 . 00000 , 8.00000 
6 . 00000 , 6 . 00000 , 8.0000 
6 . 00000 , 6 . 00000 , 12.00000 


PG: plate 2 
4, 0 

6 . 00000 , 6 . 00000 , 
6.00000, 3.00000, 
6.00000, 3.00000, 
6 . 00000 , 6 . 00000 , 


12.0000 

12.0000 

0.000000 

0.000000 


PG: plate 3 
4, 0 

6 . 00000 , 6 . 00000 , 0.000000 
6.00000, 6.00000, 4.00000 
6.00000, -6.00000, 4.00000 
6 . 00000 , - 6 . 00000 , 0.000000 


PG: plate 4 
4, 0 

6 . 00000 , - 6 . 00000 , 0.000000 
6.00000, -3.00000, 0.000000 
6 . 00000 , - 3 . 00000 , 12.0000 
6 . 00000 , - 6 . 00000 , 12.0000 

SG: Source in height 6 
0 . 000000 , 0 . 000000 , 6.000000 
0.000000, 0.000000, 90.0000, 0.000000 
- 1 , 1 . 00000 , 0.000000 
1 . 00000 , 0.000000 
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PN: Near zone cut 


9 . 00000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

3 . 00000 , 0 . 000000 , 0.000000 
0 . 000000 , 0 . 000000 , 2.00000 
7 

PL: 

T 

LY 

T 

XQ: Execute the program 

EN: End 

CE: This is the input file used for pseudo-openings 

CE: right hand ok 

FR: We will use a frequency of 2.4 GHz (ISM band) 
2.40000 

US: 

0 


PG: wall 
4, 0 

6 . 00000 , - 6 . 000000 , 12.0000 
6 . 00000 , - 6 . 000000 , 0.000000 
6 . 00000 , 6 . 0000 , 0.000000 
6 . 00000 , 6 . 0000 , 12.0000 


ZD: opening-thin material slab 
1 
1 

0.0500000, 1.01000, 0.000000, 1.00000, 0.000000 

PG: door 
4, 0 

6.00000, -3.000000, 8.0000 
6.00000, -3.000000, 4.000000 
6.00000, 3.00000, 4.000000 
6.00000, 3.0000, 8.0000 


SG: Source in height 6 
0 . 000000 , 0 . 000000 , 6.00000 
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0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

PN: Near zone cut 
9 . 00000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

3 . 00000 , 0 . 000000 , 0.000000 
0 . 000000 , 0 . 000000 , 2.00000 
7 

PL: 

T 

LY 

T 

XQ: Execute the program 

EN: End 

CE: This is the input file used for one Antenna located in 
CE: the middle of the room 

CE: Rear tank plates overlap to avoid diffraction 
CE: inside them 

CM: Dimensions of the missile room: Dim{20,10,10) Origin at 
CE: the most FORWARD and PORT point of the default view 

CM: All materials would be Perfect Electric Conductors, 

CE: unless otherwise defined. 

ZD: P.E.C. 

0 

FR: We will use a frequency of 2.4 GHz (ISM band) 

2.40000 

US: (Necessary) in Wavelengths 
0 


SG: Put the source at (10,5,2) position 

10.0000, 5.00000, 2.00000 

0.000000, 0.000000, 90.0000, 0.000000 
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- 1 , 0 . 500000 , 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT at Z=2 m. 

0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 000000 , 0 . 000000 , 2.00000 
0.0500000, 0.0500000, 0.000000 
1, 400, 200 

CE: Construct the room, according to the RIGHT-HAND RULE 

CM: Since we can not simulate the opening with an adequate 
CM: material, we define five different plates for the 
CE: Forward Level. 

CM: PG: room:(20,10,10} FORWARD level 
CM: 4, 0 

CM: 0.000000, 10.0000, 0.000000 
CM: 0.000000, 10.0000, 10.0000 
CM: 0.000000, 0.000000, 10.0000 
CE: 0.000000, 0.000000, 0.000000 

CE: Construct all the plates, but the opening itself 


PM: 

12 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0 


0.000000 

10 . 0000 , 

10 . 0000 , 

0.000000 

5.50000, 

10 . 0000 , 

10 . 0000 , 

0.000000 

5.50000, 
0.000000 

4.50000, 

4.50000, 


0.000000 

0.000000 

0.300000 

0.300000 

1.80000 

1.80000 

10.0000 

10.0000 

0.300000 

1.80000 

0.300000 

1.80000 


1 

2 

3 

4 
0 
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9 

3 
6 

5 
0 

4 

11 

12 

10 
0 

10 

6 

7 

8 
0 
0 

CM: In case we want the wall without the opening (solid), 
CE: we uncomment the following command: 


CM: PG: This is the opening in the forward bulkhead 
CM: 4, 0 

CM: 0.000000, 4.50000, 0.300000 
CM: 0.000000, 4.50000, 1.80000 
CM: 0.000000, 5.50000, 1.80000 
CE: 0.000000, 5.50000, 0.300000 

PG: room: (20,10,10) GROUND level 
4, 0 

0 . 000000 , 0 . 000000 , 0.000000 
20 . 0000 , 0 . 000000 , 0.000000 
20.0000, 10.0000, 0.000000 
0 . 000000 , 10 . 0000 , 0.000000 


PG: room:(20,10,10) UPPER level 
4, 0 

0 . 000000 , 10 . 0000 , 10.0000 
20 . 0000 , 10 . 0000 , 10.0000 
20 . 0000 , 0 . 000000 , 10.0000 
0 . 000000 , 0 . 000000 , 10.0000 


PG: room: (20,10,10) AFT level 
4, 0 

20 . 0000 , 0 . 000000 , 0.000000 
20 . 0000 , 0 . 000000 , 10.0000 
20 . 0000 , 10 . 0000 , 10.0000 
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20.0000, 10.0000, 0.000000 


PG: PORT BULKHEAD 
4, 0 

0 . 000000 , 0 . 000000 , 0.000000 
0 . 000000 , 0 . 000000 , 10.0000 
20 . 0000 , 0 . 000000 , 10.0000 
20 . 0000 , 0 . 000000 , 0.000000 

PG: STARBOARD BULKHEAD 
4, 0 

0 . 000000 , 10 . 0000 , 0.000000 
20 . 0000 , 10 . 0000 , 0.000000 
20 . 0000 , 10 . 0000 , 10.0000 
0 . 000000 , 10 . 0000 , 10.0000 

CE: First OVERHEAD: This is an intermediate floor at z 

CE:level 


PG: level:(0,0,5) first walking floor 
4, 0 

0.000000, 0.000000, 5.00000 
0 . 000000 , 10 . 0000 , 5.00000 
20 . 0000 , 10 . 0000 , 5.00000 
20 . 0000 , 0 . 000000 , 5.00000 

CE: Construct all the plates, but the opening itself 


CM: 

PG: RIGHT 


CM: 

4, 0 


CM: 

0.000000, 5.50000, 

5.00000 

CM: 

0.000000, 10.0000, 

5.00000 

CM: 

20.0000, 10.0000, 

5.00000 

CM: 

CM: 

20.0000, 5.50000, 

5.00000 

CM: 

PG: MIDDLE 


CM: 

4, 0 


CM: 

15.7500, 5.51000, 

5.00000 

CM: 

20.0000, 5.51000, 

5.00000 

CM: 

20.0000, 4.49000, 

5.00000 

CM: 

CM: 

15.7500, 4.49000, 

5.00000 

CM: 

PG: MIDDLE 


CM: 

4, 0 


CM: 

0.000000, 5.51000, 

5.00000 

CM: 

14.7500, 5.51000, 

5.00000 
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CM: 14.7500, 4.49000, 5.00000 
CM: 0.000000, 4.49000, 5.00000 
CM: 

CM: ?G: LEFT 
CM: 4, 0 

CM: 0.000000, 0.000000, 5.00000 
CM: 0.000000, 4.50000, 5.00000 
CM: 20.0000, 4.50000, 5.00000 
CE: 20.0000, 0.000000, 5.00000 

CM: In case we want the overhead without the opening 
CM:(solid), 

CE: we uncomment the following command: 

CM: PG: This is an opening in the overhead 
CM: 4, 0 

CM: 15.7500, 4.50000, 5.00000 
CM: 14.7500, 4.50000, 5.00000 
CM: 14.7500, 5.50000, 5.00000 
CE: 15.7500, 5.50000, 5.00000 


CE: These are the 6 missile tubes 

CG: First (FORWARD PORT) missile tube 
6 . 00000 , 3 . 00000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 10.0000, 90.0000 

CG: Second (FORWARD STARBOARD) missile tube 
6.00000, 7.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 10.0000, 90.0000 


RR: Make the remaining symmetric tubes. First shift 
4.00000, 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 


CG: Third (PORT) missile tube 
6.00000, 3.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 10.0000, 90.0000 
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CG: Forth (STARBOARD) missile tube 
6.00000, 7.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 10.0000, 90.0000 

RR: Make the remaining symmetric tubes. Second shift 
4.00000, 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 

CG: Fifth (PORT) missile tube 
6.00000, 3.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 10.0000, 90.0000 


CG: Sixth (STARBOARD) missile tube 
6.00000, 7.00000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
1 . 00000 , 1.00000 

0.000000, 90.0000, 10.0000, 90.0000 

CE: We have to revert to global coordinates. AGAIN 


RR: Back to global coordinates 
- 8 . 00000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 

CM: CG: Test tube at new{0,0,0). Should be placed at 
origin. CM: OK! 

CM: 0.000000, 0.000000, 0.000000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 

CM: 1.00000, 1.00000 

CE: 0.000000, 90.0000, 10.0000, 90.0000 

CE: FIRST FLOOR INTERIOR, TANKS, STAIRS, OPENING 

CM: Construct the TANKS inside the first floor. There will 
CM: be three tanks along the PORT side of the room. The 
CM: FORWARD and 

CM: AFT one will have half the height of the room (2.5m) 
and 

CM: the middle one will go up to the OVERHEAD (5m) . 

CM: Similarly 

CM: and symmetrically, will be the STARBOARD side of the 
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CM: room. 

CM: In addition, one big tank will be placed in the AFT 
part 

CE: with a height of 2.5m. 

CM: We CAN NOT consider PORT tanks, as if there was only 
one CM: large, 

CM: along the PORT side and an additional one in the center 
CM: because 

CM: of the warnings we get. For the first, 

CM: length 20m, width 0.5m and height 2.5m (half the way 
CM: up). 

CM: We only need to put the STARBOARD and UPPER plate of 
it, 

CE: since the rest are already put as walls in the room. 

PG: PORT large tank, STARBOARD-*- forward plate only 
4, 0 

0.000000, 0.500000, 0.000000 
0.000000, 0.500000, 2.50000 
5.00000, 0.500000, 2.50000 
5.00000, 0.500000, 0.000000 

PG: PORT large tank, UPPER -^forward plate only 
4, 0 

0.000000, 0.500000, 2.50000 
0.000000, 0.000000, 2.50000 
5.00000, 0.000000, 2.50000 
5.00000, 0.500000, 2.50000 

PG: PORT large tank, STARBOARD*- aft plate only 
4, 0 

15.0000, 0.500000, 0.000000 
15.0000, 0.500000, 2.50000 
20.0000, 0.500000, 2.50000 
20.0000, 0.500000, 0.000000 

PG: PORT large tank, UPPER +aft plate only 
4, 0 

15.0000, 0.500000, 2.50000 
15.0000, 0.000000, 2.50000 
20.0000, 0.000000, 2.50000 
20.0000, 0.500000, 2.50000 

PG: PORT additional tank, FORWARD plate only 
4, 0 
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5 . 00000 , 0 . 000000 , 0.000000 
5.00000, 0.000000, 5.00000 
5.00000, 1.00000, 5.00000 
5.00000, 1.00000, 0.000000 

PG: PORT additional tank, AFT plate only 
4 , 0 

15.0000, 0.000000, 0.000000 
15.0000, 1.00000, 0.000000 
15.0000, 1.00000, 5.00000 
15.0000, 0.000000, 5.00000 

PG: PORT additional tank, STARBOARD plate only 
4, 0 

5.00000, 1.00000, 0.000000 
5.00000, 1.00000, 5.00000 
15.0000, 1.00000, 5.00000 
15.0000, 1.00000, 0.000000 

PG: PORT additional tank, UPPER plate only 
4, 0 

5.00000, 0.000000, 5.00000 
15.0000, 0.000000, 5.00000 
15.0000, 1.00000, 5.00000 
5.00000, 1.00000, 5.00000 

PG: STARBOARD large tank, PORT* forward plate only 
4, 0 

0.000000, 9.50000, 0.000000 
5.00000, 9.50000, 0.000000 
5.00000, 9.50000, 2.50000 
0.000000, 9.50000, 2.50000 

PG: STARBOARD large tank, UPPER and forward plate only 
4, 0 

0.000000, 9.50000, 2.50000 
5.00000, 9.50000, 2.50000 
5.00000, 10.0000, 2.50000 
0.000000, 10.0000, 2.50000 

PG: STARBOARD large tank, PORT and aft plate only 
4, 0 

15.0000, 9.50000, 0.000000 
20.0000, 9.50000, 0.000000 
20.0000, 9.50000, 2.50000 
15.0000, 9.50000, 2.50000 
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PG: STARBOARD large tank, UPPER+ aft plate only 
4, 0 

15.0000, 9.50000, 2.50000 
20.0000, 9.50000, 2.50000 
20.0000, 10.0000, 2.50000 
15.0000, 10.0000, 2.50000 

?G: STARBOARD additional tank, FORWARD plate only 
4, 0 

5.00000, 9.00000, 0.000000 
5 . 00000 , 10 . 0000 , 0.000000 
5.00000, 10.0000, 5.00000 
5.00000, 9.00000, 5.00000 

PG: STARBOARD additional tank, AFT plate only 
4, 0 

15.0000, 9.00000, 0.000000 
15.0000, 10.0000, 0.000000 
15.0000, 10.0000, 4.99000 
15.0000, 9.00000, 4.99000 

PG: STARBOARD additional tank, PORT plate only 
4, 0 

5.00000, 9.00000, 0.000000 
15.0000, 9.00000, 0.000000 
15.0000, 9.00000, 5.00000 
5.00000, 9.00000, 5.00000 

PG: STARBOARD additional tank, UPPER plate only 
4, 0 

5.00000, 9.00000, 5.00000 
15.0000, 9.00000, 5.00000 
15.0000, 10.0000, 5.00000 
5.00000, 10.0000, 5.00000 

CM: We now construct the AFT tank, whose height is 2.5m. We 
only 

CE: need the FORWARD and the UPPER part 

PG: AFT tank FORWARD plate 
4, 0 

18.0000, 0.4500000, 0.000000 
18.0000, 0.4500000, 2.50000 
18.0000, 9.550000, 2.50000 
18.0000, 9.550000, 0.000000 
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PG: AFT tank UPPER plate 
4, 0 

18.0000, 0.4500000, 2.50000 
20.0000, 0.4500000, 2.50000 
20.0000, 9.550000, 2.50000 
18.0000, 9.550000, 2.50000 

CM: The Offset Multi-level{up and down) Ladderway will be 
put in 

CM: the center of the room. First we put a vertical plate 
at x=17 

CE: and 2=0.5. This does not count as a stair 


PG: Initial vertical plate, before first stair 
4, 0 

17.0000, 4.50000, 0.000000 
17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 
17.0000, 5.50000, 0.000000 

PG: First Stair (two plates vertical each other). 

Horizontal 

4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: First Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. First shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Second Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 
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PG: Second Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 


RR: Make the remaining symmetric stairs. Second shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 


PG: Third Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Third Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Third shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 


PG: Forth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Forth Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Forth shift 
-0.250000, 0.000000, 0.500000 
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0.000000, 0.000000, 90.0000, 0.000000 


PG: Fifth Stair {two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Fifth Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 


KR: Make the remaining symmetric stairs. Fifth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Sixth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Sixth Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Sixth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Seventh Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Seventh Stair (two vertical plates). Vertical 
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4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Seventh shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Eighth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Eighth Stair (two vertical plates). Vertical 
4,0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Eighth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Ninth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Ninth Stair, Vertical 0.1cm higher to avoid warning 

(p.161) 

4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00100 

16.7500, 5.50000, 1.00100 

CE: We have to revert to global coordinates. AGAIN 


RR: Back to global coordinates 
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2.00000, 0.000000, -4.00000 
0.000000, 0.000000, 90.0000, 0.000000 

CM: CG: Test tube at new(0,0,0). Should be placed at 
origin. CM: OK! 

CM: 0.000000, 0.000000, 0.000000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 

CM: 1.00000, 1.00000 

CS: 0.000000, 90.0000, 10.0000, 90.0000 

CM: We have already put the forward vertical plate that 
CM: surrounds 

CM: the opening. We now put the remaining three, trying to 
CE: eliminate the warnings (thus Zmax=5.Q01) 

PG: REST of the frame that surrounds the opening: LEFT 
4, 0 

14.7500, 4.50000, 5.00100 

15.7500, 4.50000, 5.00100 

15.7500, 4.50000, 4.70000 

14.7500, 4.50000, 4.70000 

PG: REST of the frame that surrounds the opening: RIGHT 
4, 0 

14.7500, 5.50000, 5.00100 

14.7500, 5.50000, 4.70000 

15.7500, 5.50000, 4.70000 

15.7500, 5.50000, 5.00100 

PG: REST of the frame that surrounds the opening: AFT 
4, 0 

15.7500, 4.50000, 5.00100 

15.7500, 5.50000, 5.00100 

15.7500, 5.50000, 4.70000 

15.7500, 4.50000, 4.70000 

CM: We should now "attach” a specific material with certain 
CM: characteristics in the overhead of the first floor in 
CM: order to simulate the opening where the stairs end! 

CM: Center 

CM: of the opening: 15 5 5 From x=14->15 y=4.5->5.5 z=5. 

CM: THIS UNORTHODOX METHOD DOES NOT WORK PROPERLY, SO WE 
CM: CHOOSE 

CM: TO WORK ANALYTICALLY, BY CONSTRUCTING EACH OTHER PART 
OF CM: THE 
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CE: OVERHEAD BUT THE OPENING ITSELF AND THEN "GLUE" THEM 
TOGETHER! 


CM: ZD: Simulated opening. Permitivity 1.01CM: 1 
CM: 1 

CM: 0.0500000, 1.01000, 0.000000, 1.00000, 0.000000 
CM: 

CM: PG: This is an opening in the overhead 
CM: 4, 0 

CM: 15.7500, 4.50000, 5.00000 
CM: 14.7500, 4.50000, 5.00000 
CM: 14.7500, 5.50000, 5.00000 
CM: 15.7500, 5.50000, 5.00000 
CM: 

CM: ZD: Back again to P.E.C. 

CE: 0 


CM: We now try to construct the HANDRAILS as cylinders with 
CM: small radius and adequate angle. USE THE Z AXIS 
CE: POSITION- ROTATE 

CG: LEFT HANDRAIL 
17.0000, 4.50000, 0.250000 
335.000, 0.000000, 0.000000, -999.000 
0 . 100000 , 0.100000 

0.000000, 45.0000, 5.00000, 45.0000 

CG: RIGHT HANDRAIL 
17.0000, 5.50000, 0.250000 
335.000, 0.000000, 0.000000, -999.000 
0 . 100000 , 0.100000 

0.000000, 45.0000, 5.00000, 45.0000 


CM: 


CM: 


END OF FIRST FLOOR- 


CE: 


PL: 

T 

XQ: 


EN: End 
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CE: This is the input file when Antenna is located in a 
CE: random place. We will present only the commands 
different from the previous file. 

CE: Put one antenna in (14,5,1). 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT at 2=1.5 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T / 

0.000000, 0.000000, 1.50000 
0.0500000, 0.0500000, 0.000000 
1, 400, 200 


EN: End 

CE: This is the input file when two Antennas are located in 
CE: the room. We will present only the co mm a n ds different 
from the previous file 

CM: Put 2 sources in Z=4m level and observation 
CM: area at Z=lm. Rear tank plates overlap to avoid 
CE: diffraction inside the tank 


SG: Put the source at (14,5,4) position 
14.0000, 5.00000, 4.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


SG: Put ANOTHER source at (6,5,4) position 
6.0000, 5.00000, 4.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
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VN: VOLUMETRIC PLOT at Z=I. 0 in. 

0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.000000, 0.000000, I.00000 
0.0500000, 0.0500000, 0.000000 
1, 400, 200 

EN: End 

CE: This is the input file for room to room case. This is 
a different file. We will present only the commands 
different from the previous file. 

ZD: P.E.C. 

0 

FR: We will use a frequency of 2.4 GHz (ISM band) 

2.40000 

US: (Necessary) in Wavelengths 
0 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

CM: Antenna is located in the first floor, but 
CM: observervation 

CM: points are in the second floor (1.5m above the floor). 
CM: Recall that the problem with the plate edges 
diffraction 
CE: still exists 

VN: VOLUMETRIC PLOT at Z=6.5 m 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

1.00000, 1.00000, 6.50000 
0.0500000, 0.0500000, 0.000000 
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1, 400, 200 


CE: Construct the ROOM, according to the RIGHT-HAND RULE 
(SAME AS BEFORE) 

CM: - 

CM; -END OF FIRST FLOOR- 

CE: - 


CE: SECOND FLOOR {MISSILES AND FRAME ALREADY PUT) 

CM: SECOND OVERHEAD: This is an Overhead at z=10 level. We 
CM: will 

CE: use the first construction, shifted by 5m along Z axis 

CM: PG: level: (0,0,10) -SECOND OVERHEAD- 

CM: Theoretically- 

CM: 4, 0 

CM: 0.000000, 0.000000, 10.00000 
CM: 0.000000, 10.0000, 10.00000 
CM: 20.0000, 10.0000, 10.00000 
CE: 20.0000, 0.000000, 10.00000 

CE: Construct all the plates, but the opening itself 


PG: RIGHT 
4, 0 

0.000000, 5.50000, 10.00000 
0 . 000000 , 10 . 0000 , 10.00000 
20 . 0000 , 10 . 0000 , 10.00000 
20.0000, 5.50000, 10.00000 

PG: MIDDLE 
4, 0 

15.7500, 5.51000, 10.00000 
20.0000, 5.51000, 10.00000 
20.0000, 4.49000, 10.00000 

15.7500, 4.49000, 10.00000 

PG: MIDDLE 
4, 0 

0.000000, 5.51000, 10.00000 

14.7500, 5.51000, 10.00000 
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14.7500, 4.49000, 10.00000 
0.000000, 4.49000, 10.00000 

PG: LEFT 
4, 0 

0 . 000000 , 0 . 000000 , 10.00000 
0.000000, 4.50000, 10.00000 
20.0000, 4.50000, 10.00000 
20 . 0000 , 0 . 000000 , 10.00000 


CM: In case we want the overhead without the 
OPENING(solid), 

CE: we uncomment the following command: 

CM: PG: This is an OPENING in the second overhead 
CM: 4, 0 

CM: 15.7500, 4.50000, 10.00000 
CM: 14.7500, 4.50000, 10.00000 
CM: 14.7500, 5.50000, 10.00000 
CE: 15.7500, 5.50000, 10.00000 

CM: The SECOND Offset Multi-level(up and down) Ladderway 
CM: will be put in 

CM: the aft part of the room. First we put a vertical plate 
CM: at x=17 

CE: and z=0.5. This does not count as a stair 

RR: Use the first construction for the SECOND LADDERWAY 

0.000000, 0.000000, 5.00000 

0.000000, 0.000000, 90.0000, 0.000000 

CM: PG: Initial vertical plate, before first stair Warnings 
CM: 4, 0 

CM: 17.0000, 4.50000, 0.000000 
CM: 17.0000, 4.50000, 0.500000 
CM: 17.0000, 5.50000, 0.500000 
CE: 17.0000, 5.50000, 0.000000 

PG: First Stair (two plates vertical each other).Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 
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PG: First Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. First shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Second Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Second Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Second shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Third Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Third Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Third shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 




PG: Forth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Forth Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Forth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 


PG: Fifth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Fifth Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Fifth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Sixth Stair (two vertical plates). Horizontal 
4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 


PG: Sixth Stair (two vertical plates). Vertical 
4, 0 


135 



16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

R R: Make the remaining symmetric stairs. Sixth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 


PG: Seventh Stair 
4, 0 

(two vertical 

plates). 

17.0000, 4.50000, 

0.500000 


17.0000, 5.50000, 

0.500000 


16.7500, 5.50000, 

0.500000 


16.7500, 4.50000, 

0.500000 


PG: Seventh Stair 
4, 0 

(two vertical 

plates). 

16.7500, 5.50000, 

0.500000 


15.7500, 4.50000, 

0.500000 


16.7500, 4.50000, 

1.00000 


16.7500, 5.50000, 

1.00000 



Horizontal 


Vertical 


RR: Make the remaining symmetric stairs. Seventh shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 

PG: Eighth Stair (two vertical plates). Horizontal 
4,0 ‘ 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Eighth Stair (two vertical plates). Vertical 
4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00000 

16.7500, 5.50000, 1.00000 

RR: Make the remaining symmetric stairs. Eighth shift 
-0.250000, 0.000000, 0.500000 
0.000000, 0.000000, 90.0000, 0.000000 


PG: Ninth Stair (two vertical plates). Horizontal 
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4, 0 

17.0000, 4.50000, 0.500000 
17.0000, 5.50000, 0.500000 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

PG: Ninth Stair, Vertical. 0.1cm higher to avoid warnings 

(p.161) 

4, 0 

16.7500, 5.50000, 0.500000 

16.7500, 4.50000, 0.500000 

16.7500, 4.50000, 1.00100 

16.7500, 5.50000, 1.00100 

CE: We have to revert to global coordinates. AGAIN 


RR: Back to global coordinates. Z TRANSLATIONS HAS CHANGED 
2 . 00000 , 0 . 000000 , - 9.00000 
0.000000, 0.000000, 90.0000, 0.000000 


CM: CG: Test tube at new(0,0,0). Should be placed at 
origin. CM: OK! 

CM: 0.000000, 0.000000, 0.000000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 

CM: 1.00000, 1.00000 

CS: 0.000000, 90.0000, 10.0000, 90.0000 


CM: We have already put the forward vertical plate that 
surrounds 

CM: the opening. We now put the remaining three, trying to 
CE: eliminate the warnings (thus Zmax=5.001) 


PG: REST of the frame that surrounds the opening: LEFT 
4, 0 

14.7500, 4.50000, 10.00100 

15.7500, 4.50000, 10.00100 

15.7500, 4.50000, 9.70000 

14.7500, 4.50000, 9.70000 


PG: REST of the frame that 
4, 0 

14.7500, 5.50000, 10.00100 

14.7500, 5.50000, 9.70000 

15.7500, 5.50000, 9.70000 

15.7500, 5.50000, 10.00100 


surrounds the opening: RIGHT 
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PG: REST of the frame that surrounds the opening: AFT 
4, 0 

15.7500, 4.50000, 10.00100 

15.7500, 5.50000, 10.00100 

15.7500, 5.50000, 9.70000 

15.7500, 4.50000, 9.70000 

CM: ZD: Simulated opening. Permitivity 1.01 
CM: 1 
CM: 1 

CM: 0.0500000, 1.01000, 0.000000, 1.00000, 0.000000 
CM: 

CM: PG: This is an opening in the overhead 
CM: 4, 0 

CM: 15.7500, 4.50000, 5.00000 
CM: 14.7500, 4.50000, 5.00000 
CM: 14.7500, 5.50000, 5.00000 
CM: 15.7500, 5.50000, 5.00000 
CM: 

CM: ZD: Back again to P.E.C. 

CE: 0 


CM: We now try to construct the HANDRAILS as cylinders with 
CM: small radius and adequate angle. USE THE Z AXIS 
CE: POSITION- ROTATE 

RR: Use the first construction for the SECOND LADDERWAY 
0.000000, 0.000000, 5.00000 
0.000000, 0.000000, 90.0000, 0.000000 

CG: LEFT HANDRAIL 
17.0000, 4.50000, 0.250000 
335.000, 0.000000, 0.000000, -999.000 
0 . 100000 , 0.100000 

0.000000, 45.0000, 5.00000, 45.0000 

CG: RIGHT HANDRAIL 
17.0000, 5.50000, 0.250000 
335.000, 0.000000, 0.000000, -999.000 
0 . 100000 , 0.100000 

0.000000, 45.0000, 5.00000, 45.0000 


CE: We have to revert to global coordinates. AGAIN 

RR: Back to global coordinates. Z GOES 5m DOWN: 
0.00000, 0.000000, -5.00000 
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0.000000, 0.000000, 90.0000, 0.000000 

CM: CG: Test tube at (0,0,0). Should be placed at origin. 
CM: OK! 

CM: 0.000000, 0.000000, 0.000000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 

CM: 1.00000, 1.00000 

CE: 0.000000, 90.0000, 10.0000, 90.0000 

LY: 

T 

CM: LP: See data in the out file the same with PL command 
CE: T 

PL: 

T 

XQ: 

EN: End 


The following, are the input sub-files for the truncation 
of the CE: room. Each room, as well as the free space 
volume, is truncated into several CELLS, in order to avoid 
observation points inside the numerous obstacles. For the 
shake of simplicity, (we would need at least 700 pages for 
all the files), we will only present the different c omm a nd s 
for antenna and observation points placements. 

ROOM DECOMPOSITION INTO CELLS (FIRST LOCATION) 

CELL 1 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at 2=1.5 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

0.100000, 1.20000, 1.50000 


139 



0.250000, 0.250000, 0.000000 
.1, 71, 3 

CELL 2 

SG: Put the source at (14,5,1 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=l. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 

T 

0.100000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 19, 7 

CELL 3 

SG: Put the source at (14,5,1 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=l. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 

T 

7.10000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 

I, 7, 7 

CELL 4 

SG: Put the source at (14,5,1) 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 

T 

II. 1000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 


position 

.000000 

m. 

.000000 


position 

.000000 

m. 

.000000 


position 

.000000 

m. 

.000000 
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1, 7, 7 

CELL 5 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at 2=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 11, 7 

CELL 6 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0/000000 
T 

0.100000, 4.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 71, 7 

CELL 7 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 19, 7 
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CELL 8 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 

I, 7, 7 

CELL 9 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

II. 1000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 10 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 11, 7 
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CELL 11 

SG: Put the source at {14,5,1} position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

0.100000, 8.20000, 1.50000 
0.250000, 0.250000, 0.000000 

1, 71, 3 

FREE SPACE 
CELL 1 


SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
I.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 1.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


CELL 2 

SG: Put the source at {14,5,1} position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 2.20000, 1.50000 
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0.250000, 0.250000, 0.000000 
1, 19, 7 


LY: 


CELL 3 

SG: Put the source at {14,5, 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at 2=1. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 

T 

7.10000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 4 


SG: Put the source at (14,5, 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=l. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 

T 

11.1000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 5 

SG: Put the source at (14,5,1) 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, C 


position 

.000000 

m. 

.000000 


position 

.000000 

m. 

.000000 


position 

.000000 
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-1, 0.500000, 0.000000 

1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

15.1000, 2.20000, 1.50000 
0.250000, 0.250000, 0.000000 

1, 11, 7 


CELL 6 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 4.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 71, 7 


EN 

CELL 7 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 
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1, 19, 7 


CELL 8 


SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

W: VOLUMETRIC PLOT 1 at Z=1.5 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

7.10000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 

1, 7, 7 


CELL 9 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11.1000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 10 


SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VM: VOLUMETRIC PLOT 1 at Z=1.5 m. 
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0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 11, 7 


CELL 11 

SG: Put the source at (14,5,1) position 
14.0000, 5.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1.5 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 8.20000, 1.50000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


ROOM DECOMPOSITION INTO CELLS (SECOND LOCATION) 
CELL 1 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1,, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2 m. 

0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 1 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


CELL 2 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
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-1, 0.500000, 0.000000 

1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1. at Z=2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 2 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 

CELL 3 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2m 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 2.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 4 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 2 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 5 

SG: Put the source at (2,2,2) position 

2 . 00000 , 2 . 00000 , 2.00000 
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0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 2.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 11, 7 

CELL 6 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2m. 

0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 4.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 71, 7 

CELL 7 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 6 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 

CELL 8 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
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-1, 0.500000, 0.000000 

1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2 in. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 6.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 9 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 6 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 10 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1 , 11 , 7 

CELL 11 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
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1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2xn 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, I 
T 

0 . 100000 , 8 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 71, 3 

FREE SPACE 
CELL 1 

SG: Put the source at (2,2,2) 


2.00000, 

2.00000, 2 

.00000 

0.000000, 

0.000000, 

90.0000, i 

-1, 0.500000, 0.000000 

1.00000, 

0.000000 


VN: VOLUMETRIC PLOT 

1 at Z=2m 

0.000000, 

0.000000, 

0.000000 

0.000000, 

T 

0.000000, 

90.0000, i 

0.100000, 

1.20000, 

2.00000 

0.250000, 
1, 71, 3 

0.250000, 

0.000000 


CELL 2 



SG: Put the source 

at (2,2,2) 

2.00000, 

2.00000, 2 

.00000 

0 . 000000 , 

0 . 000000 , 

90.0000, 

-1, 0.500000, 0.000000 

1 . 00000 , 

0.000000 


VN: VOLUMETRIC PLOT 

1 at Z=2m 

0 . 000000 , 

0 . 000000 , 

0.000000 

0 . 000000 , 

T 

0 . 000000 , 

90.0000, 

0.100000, 

2.20000, 

2.00000 

0.250000, 
1, 19, 7 

0.250000, 

0.000000 


CELL 3 


.000000 


position 

i.OOOOOO 

i.000000 


position 

i.OOOOOO 

1.000000 
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SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 2.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 4 


SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 2 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 5 


SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 
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15.1000, 2.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 11, 7 


CELL 6 

SG: Put the source at (2,2,2) 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, i 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2m 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, I 
T 

0.100000, 4.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 71, 7 


CELL 7 

SG: Put the source at (2,2,2) 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, i 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2m 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, i 
T 

0 . 100000 , 6 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 


CELL 8 

SG: Put the source at (2,2,2) 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 


position 

i.OOOOOO 

i.OOOOOO 


position 

i.OOOOOO 

i.OOOOOO 


position 

i.OOOOOO 
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1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 6.20000, 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 9 


SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 6 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 10 


SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2ffi. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 2.00000 
0.250000, 0.250000, 0.000000 
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1, 11, 7 


CELL 11 

SG: Put the source at (2,2,2) position 
2 . 00000 , 2 . 00000 , 2.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 8 . 20000 , 2.00000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


ROOM DECOMPOSITION INTO CELLS (THIRD LOCATION) 
CELL 1 

SG: Put the source at {2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2.2 m. 

0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 1 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


CELL 2 

SG: Put the source at (2,2,2.2) position 

2 . 00000 , 2 . 00000 , 2.20000 

0.000000, 0.000000, 90.0000, 0.000000 

-1, 0.500000, 0.000000 

1 . 00000 , 0.000000 
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VN: VOLUMETRIC PLOT 1 at Z=2.2 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

0 . 100000 , 2 . 20000 , 2.200000 

0.250000, 0.250000, 0.000000 

1, 19, 7 

CELL 3 

SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2.2m 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 2.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 4 

SG: Put the source at (2,2,2.2} position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2.2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 2 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 5 

SG: Put the source at {2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
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VN: VOLUMETRIC PLOT 1 at Z=2.2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 2.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 11, 7 

CELL 6 

SG: Put the source at {2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2.2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0.100000, 4.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 71, 7 

CELL 7 

SG: Put the source at {2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2.2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 6 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 19, 7 

CELL 8 

SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
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1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2.2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 6.20000, 2.200000 
0.250000, 0.250000, 0.000000 

I, 7, 7 

CELL 9 

SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2.2m. 

0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

II . 1000 , 6 . 20000 , 2.20000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 10 

SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2.2 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 11, 7 

CELL 11 

SG: Put the source at {2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 
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VN: VOLUMETRIC PLOT 1 at Z=2 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 
T 

0 . 100000 , 8 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 71, 3 

FREE SPACE 
CELL 1 


SG: Put the source at (2,2,2 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 
T 

0 . 100000 , 1 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


CELL 2 

SG: Put the source at (2,2,2 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 
T 

0 . 100000 , 2 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 19, 7 


CELL 3 


2m. 

0.000000 


2) position 
0.000000 

2m. 

0.000000 


2) position 
0.000000 

2m. 

0.000000 


159 




SG: Put the source at (2,2,2 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 
T 

7.10000, 2.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 4 


SG: Put the source at (2,2,2 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 
T 

11 . 1000 , 2 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 5 


SG: Put the source at (2,2,2 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 
T 

15.1000, 2.20000, 2.200000 


.2) position 
0.000000 

.2m. 

0.000000 


2} position 
0.000000 

2 m. 

0.000000 


2) position 
0.000000 

2m. 

0.000000 
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0.250000, 0.250000, 0.000000 
1, 11, 7 


CELL 6 

SG: Put the source at (2,2,2 


2.00000, 

2.00000, 2 

.20000 

0.000000, 

0.000000, 

90.0000, 

-1, 0.500000, 0.000000 

1.00000, 

0.000000 


VN: VOLUMETRIC PLOT 

1 at Z=2 

0.000000, 

0.000000, 

0.000000 

0.000000, 

T 

0.000000, 

90.0000, 

0.100000, 

4.20000, 

2.200000 

0.250000, 
1, 71, 7 

0.250000, 

0.000000 

CELL 7 




SG: Put the source 

at (2,2,2 

2.00000, 

2.00000, 2 

.20000 

0.000000, 

0.000000, 

90.0000, 

-1, 0.500000, 0.000000 

1.00000, 

0.000000 


VM: VOLUMETRIC PLOT 

1 at Z=2 

0.000000, 

0.000000, 

0.000000 

0.000000, 

0.000000, 

90.0000, 

T 



0.100000, 

6.20000, 

2.00000 

0.250000, 

0.250000, 

0.000000 

1, 19, 7 




CELL 8 

SG: Put the source at (2,2,2 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


2) position 
0.000000 

2m. 

0.000000 


2) position 
0.000000 

2m. 

0.000000 


.2) position 
0.000000 
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VN: VOLUMETRIC PLOT 1 at Z=2.2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7.10000, 6.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 9 


SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2.2m. 

0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 6 . 20000 , 2.200000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 10 


SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=2.2m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 2.200000 
0.250000, 0.250000, 0.000000 
1, 11, 7 
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CELL 11 


SG: Put the source at (2,2,2.2) position 
2 . 00000 , 2 . 00000 , 2.20000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=2.2 m. 

0 . 000000 , 0 . 000000 , 0.000000 

0 . 000000 , 0 . 000000 , 90 . 0000 , 0.000000 

T 

0 . 100000 , 8 . 20000 , 2.200000 

0.250000, 0.250000, 0.000000 

1, 71, 3 


EN: END 

ROOM DECOMPOSITION INTO CELLS (FORTH LOCATION) 
CELL 1 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0 . 000000 , 0 . 000000 , 90 . 0000 , 0.000000 
- 1 , 0 . 500000 , 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1 Hi. 

0 . 000000 , 0 . 000000 , 0.000000 
0 . 000000 , 0 . 000000 , 90 . 0000 , 0.000000 
T 

0 . 100000 , 1 . 20000 , 1.00000 
0 . 250000 , 0 . 250000 , 0.000000 
1, 71, 3 


CELL 2 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=1 m. 
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0 . 000000 , 0 . 000000 , 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

T 

0 . 100000 , 2 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 

CELL 3 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

7 . 10000 , 2 . 20000 , 1.00000 

0.250000, 0.250000, 0.000000 

I, 7, 7 

CELL 4 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

II . 1000 , 2 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 5 

SG: Put the source at (4,3,1) position 
4 . 00000 , 3 . 00000 , 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z-l m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
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15.1000, 2.20000, 1.00000 
0.250000, 0.250000, 0.000000 
1, 11, 7 

CELL 6 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 4 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 71, 7 

CELL 7 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=im. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 6 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 

CELL 8 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 
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7.10000, 6.20000, 1.00000 
0.250000, 0.250000, 0.000000 

I, 7, 7 

CELL 9 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

II . 1000 , 6 . 20000 , 1.0000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 10 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 1.00000 
0.250000, 0.250000, 0.000000 
I, 11, 7 

CELL 11 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 8 . 20000 , 1.00000 
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0.250000, 0.250000, 0.000000 
1 , 71, 3 

FREE SPACE 
CELL 1 


SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=im. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 1 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


CELL 2 


SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 


VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 2 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 


CELL 3 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
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1.00000, 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
•T 

7.10000, 2.20000, 1.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 4 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=1 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 2 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 5 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 2.20000, 1.00000 
0.250000, 0.250000, 0.000000 
1, 11, 7 


CELL 6 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
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1.00000, 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 

T 

0.100000, 4.20000, 1.00000 
0.250000, 0.250000, 0.000000 
1, 71, 7 


CELL 7 

SG: Put the source at (4,3,1) 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, < 
-1, 0.500000, 0.000000 
1.00000, 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, I 
T 

0 . 100000 , 6 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 19, 7 


CELL 8 

SG: Put the source at (4,3,1) 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, i 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, I 
T 

7.10000, 6.20000, 1.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 


CELL 9 

SG: Put the source at (4,3,1) 
4.00000, 3.00000, 1.00000 


.000000 


position 

i.OOOOOO 

i.OOOOOO 


position 

i.OOOOOO 

1.000000 


position 
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0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1.00000, 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

11 . 1000 , 6 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 7, 7 

CELL 10 

SG: Put the source at {4,3,1) position 
4 . 00000 , 3 . 00000 , 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 
VN: VOLUMETRIC PLOT 1 at Z=lm. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

15.1000, 6.20000, 1.00000 
0.250000, 0.250000, 0.000000 
1, 11, 7 

CELL 11 

SG: Put the source at (4,3,1) position 
4.00000, 3.00000, 1.00000 
0.000000, 0.000000, 90.0000, 0.000000 
-1, 0.500000, 0.000000 
1 . 00000 , 0.000000 

VN: VOLUMETRIC PLOT 1 at Z=1 m. 
0 . 000000 , 0 . 000000 , 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

0 . 100000 , 8 . 20000 , 1.00000 
0.250000, 0.250000, 0.000000 
1, 71, 3 


EN: END OF FORTH LOCATION 
End of all NEC-BSC files 
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APPENDIX C. MATHCAD FILES 

MathCAD file used in the comparison with free space 
(Complex Metallic Room) 


Plot of the BSC Results 
Read the Data 

NEC_BSC_data := READPRN("C:/Temp/thesis/complex_six pIate_box_test.fi> 
M := rows(NEC_BSC_data) 

M = 7.2xl0 3 
N := cols(NEC_B SC_data) 

N = 9 

phi := NEC_BSC_data ® 



NEC_BSC_data ^ 



n := 0..M - 1 



14.857 | 20.731 I 15.877 I 18.91 I 18.512 I 18.312 13.591 10.528 26.676 31.365 


N.x is the # points in the outer loop 


rows(E tot ) = 7.2 x 10 5 
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N y := 


rows! 


(Etot) 


N x 


Ny = 60 
q := 0.. Ny - 1 


startq := q-N x 
stopq := startq + N x - 1 

test ^ := submatri^Etot,startq,stop q ,0,0) 

rows(test) = 120 

cols(test) = 60 

- > 

Etot_dB : = 20-log(E toc ) 

test_dB ^ := submatri)(^ 0 t_dB, startq, stopq, 0, o) 

max(test_dB) = 39.564 

mirftest_dB) = -195.229 

max_field_d&= max(test_dB) 

dyn_range:= 23 

field_limit= max_field_dB- dyn_range 
p := 0..rows(test_dB) - 1 

fiel4 >q := if(test_dBp >q > fieldJimj£est_dBp ?q ,fieldJim^t 
test_lin ^ := submatrii( , startq, stopq, 0, o) 
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END of MathCAD file used in the comparison with free space 
(Complex Metallic Room) 

MathCAD file used in the comparison with free space (Room 
With Two Antennas) 


Plot of the BSC Results 
Read the Data 

NEC_BSC_data := READPRNrC:/Temp/thesis/fl_2ant_diff_zaxis.#> 
M := rows(NEC_B SC_data) 

M = 8x 10 4 


N 

N 


cols(NEC_BSC_data) 


phi := NEC_BSC_data 
phi T = 


< 2 » 



0 

1 

2 

3 : 

■■Mi 

'5 


7 

8 

9 

q 

1 

~T 

1 

1 

1 

1 

1 

1 

1 

1 


f NEC_BSC_data ® 


Er ; = 


10 


20 
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Riis 


c ‘a* 

4 . 


0 

0.015 

j 1.179 

2.438 

2.097 

1.952 

2.106 


Ee:= 

- 

NEC_BSC_data ^ 

10 20 

II 

yjg," 

SSI 


1 

^2 

^ 3 

4 

.'.5,, 


0.024 

0.029 

0.299 

0.606 

0.813 

1.023 


___ 


NEC_BSC_data ^ 

00 

Eq := 10 

II 



US 


SMI 

% 

tte.&Av 


§g| 

0.438 

0.219 

0.406 

0.241 

0.327 

0.201 


. \ 




E=ot ~ ^(E r ) 2 +(Ee) 2 -i-(ej ,) 2 


n := 0.. M - 1 



5^8 

1 


■ i 3 ---' 

£X1I 

ate 



' /fe' 


|&9.^ 

ii; 

0.439 

1.199 

| 2.489 

1 2.196 

2.14 

2.35 

1.115 

2.199 

3.363 

1.961 


rows(Etot) = 8 x 10 4 
rows(Etot) 


N y = 400 
q := 0.. N y - 1 
startq := q-N x 
stopq := startq + N x - 1 

test ^ := submatri^Et 0t , startq, stopq, 0,0) 

rows(test) = 200 
cols(test) = 400 

Etot_dB •= 20-Iog(E tot ) 

test_dB ^ := submatri^Eto^dB ? startq, stopq, 0, o) 

max(test_dB) = 38.777 
mir(test_dB) = -195.229 
max_field_dR= max(test_dB) 
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When dynamic range is 30 db, we get some errors, thus we 
set the following: 
dyn_range := 25 

field_limit= max_field_dB- dyn_range 
p := 0..rows(test_dB) - 1 

field,, q := if(test_dBp, q > fieldJimitest_dBp, q ,field_lim)t 
test_lin ^ := submatri)( Et 0 t,starlq, stopq, 0, o) 
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END of Ma tihCAD file used in the comparison with free space 
(Room With Two Antennas) 

MathCAD file used in the comparison with free space 
(Symmetric placement of Antenna) 

Plot of the BSC Results 
Read the Data 

islECURSCJlatl 
M := rows(NEC_BSC_data) 

M = SxlO 4 

N := cols(NEC_BSC_data) 

N = 9 

phi := NECJBSC_data ® 
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Ny = 400 
q := O..Ny - 1 
startq := q-N x 
stopq := starlq + N X - 1 

test ^ := submatri^Etot,start^,stop q ,O.O) 

rows(test) = 200 
cols(test) = 400 
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Etot_dB := 20-log(& ot ) 

test_dB ^ := submatri^Etot_dB,stariq.stop q ,0.0) 

max(test_dB) = 65.84 

miiftest_dB) = -195.229 

maxjfield_dR= max(test_dB) 

dyn_range:= 41 

fieldjimit= max_field_dB- dyn_range 
p := 0.. rows(test_dB) - 1 

fiel^ >q := if(test_dBp q > field_limittest_dBp, q .fieldJim)t 
test Jin ^ := submatri^E^t. stariq, stop q , 0 , o) 



[field 


END of MathCAD file -used in the comparison with free space 
(Symmetric placement of Antenna) 
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MathCAD file used in the comparison with free space (Random 
placement of Antenna) 


Plot of the BSC Results 
Read the Data 

NECJBSCLdata := RE.4DPRN("C:/remp/thesis/New_lf_r_op_an_pl_l anti* 
M := rows(NEC_BSC_data) 

M = 8 x 10 4 

N := cols(NEC_BSC_data) 

N = 9 


phi := NECJBSC_data ® 



NEC_BSC_data ^ 


E* := 10 

vJ = \ 1 <> 1 1 I g [ 3 | 4 I5 

^ 0 0.39 0.247 0.252 0.257 0.263 0.268 


E«“ 1 /( E f) 2+ ( Ee ) 2+ ( E *) : 

n:=0..M- 1 


K.x is the # points in the outer loop 
rows(Et 0t ) = 8 x 10 4 


7 

8 

9 : 

0.281 

0.286 

0.291 











rows(Etot) 

JNy I — 

N x 

N y = 400 
q := 0.. N y — 1 
startq := q-N x 
stopq := startq + N x - 1 

test ^ := submatri^Etot,startq,stopq,0,0) 

rows(test) = 200 

cols(test) = 400 

- » 

Etot_dB := 20-log(E tot ) 

test_dB ^ := submatri)(Etot_dB»startq, stopq, 0,0) 

max(test_dB) = 44.709 
min(test_dB) = -195.229 
max_field_dB= max(test_dB) 

When dynamic range is 30 db, we get some errors, thus we 
set the following: 

dyn_range := 33 

field_iimit= max_field_dB- dyn_range 
p := 0.. rows(test_dB) - 1 

fielt^q := if(test_dBp >q > fIeld_limjltest_dBp > q, fieldjim^t 
testjin ^ := submatri^Et 0t , startq, stopq, 0,0) 
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ill 


END of MathCAD file used in the 
(Random placement of Antenna) 


is: 




Sf 5F 3 




comparison with free space 
















MathCAD file used in the comparison with free space (Room 
to room scenario) 

Plot of the BSC Results 
Read the Data 

NEC_BSC_data := READPRN ("C:/Temp/thesis/New_2F_r_op_anaI_pIates.txt" 

M := rows(NEC_BSC_data) 

M = 8xl0 4 

N := coIs(NEC_BSC_data) 

N = 9 

phi := NEC_BSC_data ® 


6.51 6.5 


NEC BSC data 




MO -10 | MO -10 | MO -10 | MO -10|MO -10|MO -10 


NEC BSC data 


MO -10 |MO -10 | MO -10 | M0 -10|MO -10JMO *10 


NEC BSC data 


1-10 ' 10 

Jjk 

O 

i 

o 

1-10 - 10 

MO - 1 ° 

M 0 -10 

o 

O 


Etot :=V(Er) 2 ^(Ee)^(Eo) 
n:=0..M-l 


Six is the # points in the outer loop 
rows(Etot) = 8 x 10 4 



V?,' 

1.732-10 *10 

1.732-10 *10 










































rowset) 

N x 


Ny = 400 

q := O..Ny - 1 


startq := q-N x 

stopq := startq + N x - 1 

test ^ := submatii}(Et 0 t, startq, stopq, 0,0) 

rows(test) = 200 

cols(test) = 400 

- > 

Etot.dB •= 20-log(E t ot) 

test_dB ^ := submatri^Etot.dB.startq,stopq,0.0) 

max(test_dB) = 23.308 


min(test_dB) = -195.229 
max_fieId_dR= max(test_dB) 

When dynamic range is 30 db, we get some errors, thus we 
set the following: 

dyn_range := 75 

fieid_limit= max_field_dB- dyn_range 


p := 0..rows(test_dB) - 1 

fielc^ q := if(test_dBp, q > field_limitest_dBp jq ,field_lim)t 
testjin ^ := submatri^Et 0t , startq, stop q , 0 , o) 
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END of MathCAD file used in the comparison with free space 
(Room to room scenario) . 


MathCAD file used in calculating the excess gain/loss 
random variables (First location) 


NEC_BSC_data := READPR^dataroom" > 

NEC_BSC_data ® 

>’coord := NEC_BSC_data (l) 

NEC_BSC_data ® 

antenna , obser.points NOT in the same level 
. . .> 

P ~|w - l 4 - 0 ) 2 + (y'coord - 5.0) 2 + (wn - L5 ) 2 


free_space_data := RE ADPRN "datafiree'' ) 


f NEC BSC_data ® ^ 


4- 


10 


20 


NEC_BSC_data 


<s> 


%:= 10 


20 


NEC_BSC_d2T3 


(i> 


E,:= 10 


20 
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rows(Et 0t ) = 1.539x I0 J 


/ Q \ 

free_space_data" 1 


^cJxee-~“ 


10 


20 


free_spacc_data 


<5> 


£y_free : *~ ^ 


20 


free_space„dat 2 


<?) 


£z_free : ““ ^ 


20 


Qo:Jrce : ~ \JfexJrd~ * 0%'jfrc^~ + (^zjre<)~ 

Since the NEC-BSC dynamic range is only 30dB, we will have 
to account for this fact by "filtering 15 out all values 
which are out of that range. This is done in the following 
fashion: 


Eto idB == (201og (Ept)) 
jj:=0..rows(E touiB )-l 

lowest_allowed_Etot_value := ma^E^^g) - dyn_range 
lowest_allowed__Etot„value = 12.948 


%ot„free_dB (20‘lOg(E to t_fre^) 

lowe$uallowed_Etot_fr_val := ma ^E to tjice_dB) ” dvn_range 


lowest_aiIowed_Etot_fr_val = 7.82 

Ejotai 4B.. : =i^Etot_dB. ) < lowest_alIowed_Etot„vaiue ,lowest_allowedJEtot_value ,15^® J 

Etotai free dB *= #feot free dB < !owe$t_allowed3tot^val, lowest„ailowed3tot„fr^val ^ oxjTCt _ m 
“ " ii \ “ " ii 

Converting to linear: 

-» 

Etotal-10 20 

-> 

fr*x>uu_ftxc_d£ 

10 ^ 

Let us now create the two data matrices, namely the 
"datamatrix" which consists of the room's field strength 
values for every observed point and the ”freespacedataOl” 
which includes the corresponding data for the free space 
case. 

datamatrix:= augment (p . E^i ) 
freespacedataO 1 := augment (p, 
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At this point, owing to the nonuniform layout of the 
distance (from the antenna) points in the □ array, we will 
have to distinguish them into several different groups 
using the following approach: 
lower := floor{min(p )) 
upper := ceil(ma^p)) 

Maximum anri m-i n-i rmiTn distance for this antenna location: 
upper =15 
lower = 0 

Groups separation for this calculation (in cm): 

separation :=50 

"Step factor" for this calculation: 

./ 100 

factor := ceil - 

^ separation 

factor = 2 

Number of distance bins for this calculation: 

bins := upper-factor 

bins = 30 

, upper - lower 

h :=- 

bins 

h = 0.5 
n := 0.. bins 
k n := lower -f n*h 

It should be noted that the number of bins may be chosen in 
a fashion which best suits our desirable “distance steps”. 
For instance, should we choose to have groups points 
separated by 50cm, where the maximum distance from the 
antenna is 49m, our intuitive selection should be 98. 

We next introduce a new function which shall be utilized 
several times in our subsequent calculations. As input, it 
expects a two column matrix and an array. It looks for 
those values of the first column (i.e., in our problem, the 
distance array which fall within our made-up groups and 
then creates a new matrix that contains all those second 
column's values whose respective first column value satisfy 
our requirement. It finally returns this new matrix which, 
it should be underlined, normally has a different number of 
elements in each row. 
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build(M,bmvector) := 


Af-M^ 

B <- M ^ 
k <— binvector 
for 1 € 0.. rows(k) — 2 
Dj ^— 0 

for i€ O..rows(k) - 2 
for js 0..rows(M) - 1 
D; «- augment(Bj,Dj) if Iq < Aj < kj+j 
for i € 0.. rows(D) - 1 

D; <- submatrix(D i ,0,rows(D i ) - l,0,cols(D;) - 2) if D. * 0 
D 

It is high time we used this helping function. 
datamatrix02:= build(datamatrix, k) 
freespacedata02 

datamatrixoJ = 

The following new function is, seemingly, of high 
importance since it is the one which calculates the 
"normalized difference" random variables. It receives two 
matrices of the same dimensions and returns a new matrix 
which in each row contains the values of the random 
variable pertinent to that row's "distance group.” 


:= build(freespacedata01 ,k) 



0 

1 

2 

3 


; : 5 ^ 

6 

0 

[1,13] 

[1,35] 

[1,31] 

[1,40] 

[1,67] 

[1,93] 

[1,127] 
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difference^, Y) := 


for k e 0.. rows(X) - 1 
Ck <- 0 

f for j € O-.rows(X) - 1 
n <— cols(Xj) 

B Yj ' 
ie -0 

for is 0..n - 1 




f A„ ^ 


diff; 4 - 20-log 


A 0, 

®0,i 


«ol$(Xj) 


XI 

Cj <- stack (Cj, diff) , 

Cj < -{Cj) T 
k <- rows(Cj) 
ll— cols(Cj) 

Cj <- submatrix(Cj,0,k - 1,1,1 - l) 


The "room gain" is defined here in dB 

We now get, rid of the redundant first row zeros and then 
calculate the normalized difference variables. 
new_datamatrix02 := submatrix(datarnatrix021, rows(datamatrix02) — 2,0. cols(datamatrix02> - 1) 
new_freespacedata02 := submamx(freespacedata02 , Lrows(freespacedata02) - 2,0, cols(freespacedata02) - 1) 


norm. 


diff 1 = 


ilill 

gpgpg 

1 

2 

13lll 

l! ^ 1 


y> 6 

0 

[1,35] 

[1,35] 

[1,40] 

[1,67] 

[1,93] j 

[1,127] 

[1,134] 


form a bug for it returns an erroneous matrix. The cause of 
this weird occurrence remains unknown, however we shall 
attempt to effectively overcome this mishap by applying an 
easy trick, 
i := 0.. rows(new_datamatrix02) - 1 
rj := cols (new_datamatrix02 i ) 

corr_norm_diffj := submatrix(norm_diff i ,0,rows{norm_diff i ) - 1,0,r; - l) 


corr_norm_diff 



i|f ; 1 

|*2 

; ^ 3' 

1 '4 • 

HHi 

|§§if ^ 

0 

[1,35] [1,31] 

[1,40] 

[1,67] 

[1,93] 

[1,127] 

[1,134] 


This last matrix named "corr_norm_diff” is the one we 
intended to find. In each row of unequal element number, it 
contains a different random variable as function of our 
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"devised" groups. Let us take a closer look in one of them, 
say the 5th. 



mean(test_l 0 ) = 0.592 

THESE ARE % DIFFERENCE 

Defining the number of the bins: 

Npoims : =200 

il •” 0.. 1 

j := 0.. rows(corr_norm_diff) - 1 

upper j := ma^corr_norm_diff j) 

lowerj := min(corr_norm_diff j) 

upper_max := nuu(upper) 

lower_min := min( lower) 

upperjcnax - lower_min 
step :=--- 

‘'points 

upper_max = 13.584 
lower_min = -17.064 


step =0.153 


genjntervals „ := lower_min + step -ii 



rows(gen_intervals ) = 200 


buildUL) := 


•4 


for q € l..rows(L) - 1 
A Lj 

4ew stack (l^, A) 


Lpdf := buildl^L) 
rows(4 d f) =28 
cols(lp<ff) = 199 
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qq “CLNpdft- 1 

avg_pdf qq := meaa(L p< j f ^ ) 


^avg_pdf = 0.998 


avg pdf 


V- 




IilftlvlS 

4 

.;i^yajrtTT^ffr^. L-trtn * t. j * j ■ 

ill 

3.84-10 - 4 

0 

0 

0 

3.436-10 - 4 | 

0 
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•RT^ n of MathCAD file used in calculating the excess 
gain/loss random variables (First location) 


MathCAD file used in calculating the excess gain/loss 
random variables (Second location) 

Calculating the normalized difference random variables 
Random position STACK 

In this program we are aiming at calculating all the random 
variables of the normalized difference between the room's 
and the free space field strength (Etot and Etot_free 
respectively).We start by reading the associated data files 
generated by NEC_BSC. 

NEC_BSC_dataI := READPRKX^temp/Thesis/xla.txt") 

Ml := rows(NECJBSC_datal) 

Ml = 213 

N1 := cols (NEC_BSC_datal) 

N1 = 9 

NEC_BSC_data2 := I^ADPRM , ’C^n^^esiS/j(^^t|) 

M2 := ro ws (NECJB SC_data2) 
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M2 = 133 

N2 := cols (NEC_BSC_data2) 

N2 = 9 

NEC_BSC_data3:= READPRK"C:/temp/Thesis/x3a.txt") 
M3 := rows (NECJBSC_data3) 

M3 =49 

N3 := cols(NEC_BSC_data3) 

N3 = 9 

NEC_BSC_daia4:= READPR>TC:/temp/Thesis/x4a.txt") 
M4 := rows(NECJBSC_data4) 

M4 = 49 

N4 := cols(NEC_BSC_dara4) 

N4 = 9 

NEC_BSC_data5 := PEADPR^TC:/temp/Thesis/x5a.txt" ) 
M5 := rows (NEC_B SC_data5) 

M5 = 77 

N5 := cols(NEC_BSC_data5i) 

N5 = 9 

NEC_BSC_data6:= READPRN''C:/temp/Thesis/x6a.txf) 
M6 := rows (NEC_B SC_dataQ 
M6 = 497 

N6 := cols (NEC_BSC_data6) 

N6 = 9 

NEC_BSC_data7 := READPR^"C:/temp/Thesis/x7a.txf ) 
M7 :=rows(NEC_BSC_data7) 

M7 = 133 

N7 := cols (NEC_BSC_data7) 

N7 = 9 

NEC_BSC_data8 := READPRK"C:/ffimp/Thesis/x8a.txt") 
M8 :=rows(NEC_BSC_data8) 

M8 = 49 

N8 := cols(NEC_BSC_data8) 

N8 = 9 

NEC_BSC_data9 := READPRI^ ”C:/temp/ThesIs/x9a.txt" ) 
M9 :=rows(NEC_BSC_data9) 

M9 = 49 

N9 := cols(NEC_BSC_data9) 

N9 = 9 

NEC JBSC_datalO:= RE ADPRN( "C:/temp/Thesis/x 1 Oa.txt” ) 
M10 := rows(NEC_BSC_datalQ 
M10 = 77 

N10:= cols(NEC_BSC_datalQ 
N10 = 9 

NEC_BSC_datal 1 := READPRNCCr/temp/Thesis/xl la.txt") 
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Mil := rows(NEC_BSC_datal 1) 

Mil = 213 

Nil := cols ( NEC_BSC_datal 1) 

Nil = 9 

Xcoord- NEC_BSC_data ^ 

y coord :=NEC_BSC_data {l) 

Zcoord- NECJBSC.data <2 > 

. .— . .. > 

P ^-Ji^coord ~ -} + (ycoord ~ 2j“ f (z^ord ~ 
free_$pace_datal :=READPRN"C:/temp/Thesis/fsla.txt") 
MM1 := rows(free_space_datal) 

MM1 = 213 

NN1 := cols(free_space_datal) 

NN1 = 9 

free_space_data2 := READPRN"C:/temp/Thesis/fs2a.£xt") 
MM2 := rows(free_space_data2) 

MM2 = 133 

NN2 := cols(free_space_data2) 

NN2= 9 - 

free_space_data3 := READPRN"C^temp/Thesis/fs3a.txt" ) 
MM3 := rows(firee_space_data3 ) 

MM3 =49 

NN3:= cols(free_$pace_data3 ) 

NN3 = 9 

free_space_data4 :=READPRN"C:/temp/Thesis/fs4a.txt" ) 
MM4 := rows(free_space_data4 ) 

MM4 = 49 

NN4:= cols (free_space_data4) 

NN4= 9 

free_space_data5 := READPRN”C:/temp/Thesis/fs5a.txt'' ) 
MM5 := rows (free_space_data5 ) 

MM5 = 77 

NN5 := cols(free_space_data5) 

NN5 = 9 

free_space_data6 := READPRN "C:/temp/Thesis/fs6a.txt” ) 
MM6 := rows (free_space_data6) 

MM6 = 497 

NN6:= cols (free_space_data6) 

NN6 = 9 

free_space_data7 := READPRlN"C:/temp/Thesis/fs7a.txt” ) 
MM7 := rows (free_space_data7 ) 

MM7 = 133 



NN7 := cols(free_space_data7) 

NN7 = 9 

free_space_data8 := RE.AJDPR><"C:/temp/Thesis/fs8a.txt" ) 
MM8 := rows (free_space_data8 ) 

MM8 = 49 

NN8 := cols(free_space_data8) 

NN8 = 9 

free_$pace_data9 := READPRK T:/temp/Thesis/fs9a.t\t" ) 
MM9 := rows(free_space_data9 ) 

MM9 = 49 

NN9 := cols(free_space_data9) 

NN9 = 9 

free_space_datalO := READPRI^"C:/temp/Thesis/fslOa.txt" ) 
MM10 := rows(free_space_datalO) 

MMIO = 77 

NN10:= coIs(free_space_datalO) 

NN10= 9 

free_space_datal 1 := READPRN"C:/temp/Thesi$/f$l la.txt" ) 
MM11 :=rows(free_space_datall) 

MM11 = 213 

NN11 := cols{free_space_datal 1) 

NN11 = 9 


NECJBSC_data 


< 3 > \ 


Ex := l 10 


20 


NEC_BSC_daa ^ 


Ey := 10 


20 


NEC_BSC_daa 


<7) 


Ej- 10 


20 


Sot :=J(E x y+(E y ) 2 +(E z ) 2 
rows(E tot ) * 1.539x 1 (f 


f freejspace^data^ ^ 




10 


20 


free„space_data 


< 5 > 




20 


frec_$pace_data 


.<?> 




20 
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^tor_free - ^x_fre^ ~ ' r "** 

Since the NEC-BSC dynamic range is only 30dB, we will have 
to account for this fact by "filtering” out all values 
which are out of that range. This is done in the following 
fashion: 


E ot _ dB :=(201og(E lo; )) 



j] := 0.. rows(E tot _ dB ) - 1 

lowest_a!lowed_Etot_vaIue := ma^E^as) - dyn_range 


lowest_allowedJEtot_value = 22.107 
-. 

Eto:_free_dB := (“^l°§(^tot_fre^) 

lowest_allowed_Etot_fr_val := - dyn_range 


lowest_allowed_Etot_fir_val = 22.47 

EtoaLdB “i^jEtot dB..J < lowest_a!lowedJEtotjvaiue ,lowest_allowedJEtot_value 

Etotai_free_dB.. “ if^oUree.dB,. < lowest_ailowedJEtotJfr_val, lowest_allowed_Etot_fr_val ,Etotjree_dB ;; j 

Converting to linear: 

- > 

Etotal := 10 20 

- > , 

Etc aLfee- 10 “ 

Let us now create the two datamatrices, namely the 
"datamatrix" which consists of the room's field strength 
values for every observed point and the ”freespacedataOl" 
which includes the corresponding data for the free space 
case. 

datamatrix := augment(p jE^) 
freespacedataOl := augment(p, Etetaij&e^ 

At this point, owing to the nonuniform layout of the 
distance (from the antenna) points in the 0 array, we will 
have to distinguish them into several different groups 
using the following approach: 
lower := floor{mm(p )) 
upper := ceil(ma^p)) 

Max imum and mi nimum distance for this antenna location: 


upper = 17 
lower = 0 

Groups separation for this calculation (in cm): 

separation := 50 

"Step factor" for this calculation: 
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./ 100 

factor := ceil - 

^ separation 

factor = 2 

Number of distance bins for this calculation: 

bins := upper-factor 
bins =34 

, upper - lower 

h :=- 

bins 

h = 0.5 
n := 0.. bins 
k a := lower -f n-h 

It should be noted that the number of bins may be chosen in 
a fashion which best suits our desirable “distance steps". 
For instance, should we choose to have groups points 
separated by 50 cm, where the maximum distance from the 
antenna is 49m, our intuitive selection should be 98. 

We next introduce a new function which shall be utilized 
several times in our subsequent calculations. As input, it 
expects a two column matrix and an array. It looks for 
those values of the first column (i.e., in our problem, the 
distance array □..which fall within our made-up groups and 
then creates a new matrix that contains all those second 
column's values whose respective first column value satisfy 
our requirement. It finally returns this new matrix which, 
it should be underlined, normally has a different number of 
elements in each row. 

build(M,binvector) := A <— M ^ 

BeM^ 
k <— binvector 
for le 0..rows(k)-2 
Dj <- 0 

for i e 0.. rows® - 2 
for je 0..rows(M) - 1 
D. <=— augment (Bj,D ; ) if k. < Aj < k._j 
for i 6 0.. rows(D) - 1 

Dj <— submatrix(Di,0,rows(D i ) - 1,0, cols (D*) - 2) if D ; # 0 

D 

It is high time we used this helping function. 
datamatrix02:= build(datamatrix, k) 
freespacedata02 

datamatrixO^ = 
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:= build(freespacedata01 ,k) 






The following new function is, seemingly, of high 
importance since it is the one which calculates the 
"normalized difference" random variables. It receives two 
matrices of the same dimensions and returns a new matrix 
which in each row contains the values of the random 
variable pertinent to that row's “distance group." 


difference( X, Y) := 


for k€ 0..fows(X) - 1 

f for j € 0..rows(X) - 1 
n 4- cols(Xj) 

A 4— Xj 

Be- Yj 
I 

'for ie : 0..ri T- l 
diff 4-20-lod 


r*.-* 


m 


B, 


■OJ: 


n 4- coIs(Xj) 

Cj 4- stack(Cj.diff) 

c j 4- ( C j) T 

k <— rows(Cj) 

1 <- cols(Cj) 


Cj <- submatrix(Cj,0,k -1.1,1- l) 


The "room gain" is defined here in dB 

We now get rid of the redundant first row zeros and then 
calculate the normalized difference variables. 
new_datamatrix02:= submatrix(datamatrix02 l,rows(datamatrix03 — 2,0, cols (datamatrix02> — 1) 
new_freespacedata02 := submatrix(freespacedata02,1, rows (freespacedata02) - 2,0. cols (freespacedata02) - 1) 
norm_diff := difference( new_datamatri x02, new_freespacedata02) 


norm_( 


diff T : 


rv -,v 

0 ' 



4 

;§j 

6 

H 

[1,32] 

[1,46] [1,51] 

[1,51] 

[1,51] 

[1,51]1 

[1,51] 


Surprisingly and unexpectedly the Mathcad appears to suffer 
form a bug for it returns an erroneous matrix. The cause of 
this weird occurrence remains unknown, however we shall 
attempt to effectively overcome this mishap by applying an 
easy trick. 

i := 0.. rows(new_datamatrix02) - 1 
r ; := cols (new_datamatrix02j) 

corr_norm_<iiffi := submatrix(noniL.diff i ,0,rows{norm_<iiff i ) - 1,0,r, - l) 
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corr_norm_diff = --- ■ v --^^ . -„i .,----- 

0 I [1,32] | [1,46] 1 [1,51 ] 1 [1,47] [ [1,50] | [1,48] | [1,48] 

This last matrix named "corr_norm_diff “ xs the one we 
intended to find. In each row of unequal element number, it 
contains a different random variable as function of our 
"devised" groups. Let us take a closer look in one of them, 
say the 5th. 


0 


2 

3 4 

5 ■ 

6 

7 J 

8 

0 2.395 

-0.754 

2.494 - 

1.278 -0.751 

-0.465 

0.245 

-2.726 

-0.447 


test_0 := vcorr_nonn_diff 
- 

test_0 0 = -___—- 

ITT 2.3951 - 

mean(test_0o) = 0.338 



; ; . 2 

3 

4 

5 

6 

• 7 

8 

1.31 

-0.076 

2.766 

3.854 

-6.99 

2.108 

2.363 

5.748 


test_l := \corr_norm_diff j ___ 

, |0ll I 2 ’ I 3 I • 
oj 2.381 1.311-0.0761 2.7661 3.! 

mean(test_l 0 ) = 1.121 

THESE ARE % DIFFERENCE 

Defining the number of the bins: 

^points := 

ii . = 0.. 1 

j := 0.. rows(corr_noim_difD - 1 

upper j := xna^corr_norm_diff j) 

lowen := min(corT_norm_diff j) 

upper_max := nm( upper) 

lower_min := min( lower) 

upper max - lower_min 

step := -=_- 

^points 

upperjtnax = 13.48 
lower min = — 11.543 


step = 0.125 

gen Jntervals & := !ower_min -f step -ii 

hist(genJntervals ,conr^orao_diffj) 

jh&Z— TTT . ._' V •: 

; 3 cols(corr_norm_<iiffj) ' 

rows(genJntervals ) = 200 
buildL(L) := L, cw 4- ij 

for q e 1.. rows(L) - 1 

A Lq 1 " 

Laew ^ stack (Lj, ew , a) 

Lpdf := buildL(L) 
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qq ••= o.. -1 

avg_pdf qq := mean(L pdf ^ ) 


^avg_pdf = I 


200 







































































































































































avs^pdr =_I--- 

0 { 5 . 894-10 - 3 | 0.4741 5 .^ 

^pdfs~ ‘ 

mean_rv := ^ genjntervals qq -avg_pdf qq 
qq =0 

iroeaajrv = 1.435 
step =0.125 
avg„pdf §9 := 0.005 

^avg__pdf =0.531 

„ avs_pdf 
avgjxlf :=--- 

^avg_pdf 
^Tavg_pdf = 1 

mean_sq_rv := ^ (gen_intervals qq ) 2 -avg_pdf, 
qq=0 

mean_sq_rv = 25.604 

*> 

var„rv := mean_sqjrv - meaiL_rv~ 
varjrv = 23.545 
mean_rv = 1.435 

Defining Gaussian pdf; 

- (g<mjmervalsa- mean„rv) 2 
1 2'var rv 

gauss a := — . -e " : 


^2'71'Vdijy 

room_gain dB := genjntervals 


91 [ 92 v' 93 , 

5.734-10 -3 5.304-10 -3 7.015-10 ~ 3 
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Antenna at location 2 
























END of MathCAD file used in calculating the excess 
gain/loss random variables (Second location) 


MathCAD file used in calculating the excess gain/loss 
random variables (Third location) 


NEC_BSC_data:=READPRH"dataroom_b” ) 

^coord : = NEC_BSC_data <°> 
ycoord := NEC_BSC_data ^ 

^oord“ NEC_BSC_data ® 

antenna , obser.points in same 
—.- '. —± 

P * = -^(^coord — “)” + (y’coord — -)” + (^coord -•-) 

free_space_data := READPRI^"datafree_b” ) 


f NECJBSCJta ® ^ 


E> ; := 


10 


20 


NEC_BSC_data 


< 5 > 


E*:= 10 


20 


NEC_BSC_data 


<7> 


E*:= 10 


20 


level 
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rows(E tot ) = 1.539x 10 5 


. A .0> \ 

free.$pace_da& ' 


%cjree : ~ 


10 


20 


free„space„da& 


< 5 ) 


£yjxee : “ ^ 


20 


jxec_space_dan 


< 7 > 


Ez.free-10 


20 


%ot_free yj(^x_frc^“ * * (^Jra) 

Since the NEC-BSC dynamic range is only 30dB, we will have 
to account for this fact by "filtering" out all values 
which are out of that range. This is done in the following 
fashion: 

Etot_dB := (20-log(E tot )j 


jj:=0..rows(E tot _ dB )-1 

lowest_aUowed_Etot_value := ma?(E 0 ._ dB ) - dyn_range 
lowest_aUowedJEtot_value =22.115 

-4 


Etot_free_dB(20- log (Ejotjxe^ j 

lowest_aUowedJEtot_fir_val := ma>(E tot _ free _ dB ) - dyn_range 


Iowest_allowed_Etot_fr_va] = 22.47 

E 0t2 ] (jg rsifjjifc* < lowest_allowed_Etot_value , lowest_aIlowedJEtot_value 

Stoai free dB •- iffeot free <3B < lowest_allowed_Etot_fr_val r lowest_allowed_Etot_fr_val ,Etot &ee dB ) 

— iS \ — ~ ii ~ ~ til 

Converting to linear: 

-» 

totaI_d£ 

Eotai “ 10 20 

-» 

EtotaLfree- 10 20 

Let us now create the two datamatrices, namely the 
"datamatrix" which consists of the room's field strength 
values for every observed point and the "freespacedataOl” 
which includes the corresponding data for the free space 
case. 

datamatrix := augment (p 
freespacedaraOl := augment(p ,E rotaI _&^ 
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At this point, owing to the nonuniform layout of the 
distance (from the antenna) points in the □ array, we will 
have to distinguish them into several different groups 
using the following approach: 
lower := flooi{niin(p)) 
upper :=ceil(ma>{p)) 

Mavimrim and minimum distance for this antenna location: 


upper = 17 
lower = 0 


Groups separation for this calculation 

separation :=50 

"Step factor" for this calculation: 

( 100 'i 


factor := cei 


separation 


(in cm): 


factor = 2 

Number of distance bins for this calculation: 

bins := upper-factor 


bins = 34 


, upper - lower 

h :=- 

bins 

h = 0.5 
n := 0.. bins 


k n := lower + n-h 

It should be noted that the number of bins may be chosen in 
a fashion which best suits our desirable “distance steps". 
For instance, should we choose to have groups points 
separated by 50 cm, where the maximum distance from the 
antenna is 49m, our intuitive selection should be 98. 

We next introduce a new function which shall be utilized 
several times in our subsequent calculations. As input, it 
expects a two column matrix and an array. It looks for 
those values of the first column (i.e., in our problem, the 
distance array □..which fall within our made-up groups and 
then creates a new matrix that contains all those second 
column's values whose respective first column value satisfy 
our requirement. It finally returns this new matrix which, 
it should be underlined, normally has a different, number of 
elements in each row. 


205 



build(M, binvector ) := 


A 4— M ^ 

B 4- M ^ 
k 4 — binvector 
for le O..rows(k) -2 
Dj 4— 0 

for i € 0.. rows(k) - 2 
for j e 0..rows(M) - I 
D; 4— augment(B r D;) if k < A: < k_i 
for i e 0.. rows(D) - 1 

Dj, 4- submatrix(D i ,0,rows(D;) - 1,0, cols (Dj) - 2) if D ; * 0 
D 

It is high time we used this helping function. 
datamatrixQ2:= build(datamatrix, k) 
freespacedata02 

datamatrixO^ = 

The following new function is, seemingly, of high 
importance since it is the one which calculates the 
"normalized difference" random variables. It receives two 
matrices of the same dimensions and returns a new matrix 
which in each row contains the values of the random 
variable pertinent to that row's "distance group." 


:= buiM( freespacedataO 1. k) 



0 




■ 


6 

H 

[1,10] 

[1,32] 

[1,46] 

[1,51] 

[1,47] 

[1,50] 

[1,48] 
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difference( X,Y) := 


for ke 0..rows(X) - I 
C};<-0 

( for j € 0..rows(X) - 1 
n <r- cols(Xj) 


r\ t :— 

B «- Yj 

i 4 — 0 

for is 0.. n - 1 


diffj <— 20-log 




%.i 


n 4 — cols(Xj) 

Cj <— stack(Cj.diff) 

Cj<-(Cj) T 
k <r~ rows(Cj) 

1 ^— cols (Cj) 

Cj <r- submatrix(Cj ? 0 ? k - 1,1,1—1) 


The "room gain" is defined here in dB 

We now get rid of the redundant first row zeros and then 
calculate the normalized difference variables. 
new„datamatrix02:= submatrix(datamatrix02 l ? rows(datamatrix03 - 2,0,co!s(datamatrix0?) - 1) 
newjree$pacedata02 := submatrix(freespacedata02 , l,rows(freespacedata02) - 2 ? 0, co!s(free$pacedata02) - 1) 


norm. 


diff 1 ' = 



0 

1.- 

2 

3 

4 

Icll 

6 

0 

[1,32] 

> . * 

CO 

V— 

[1,51] 

[1,51] 

[1,51] 

[1,51] 

[1,51] 


form a bug for it returns an erroneous matrix. The cause of 
this weird occurrence remains unknown, however we shall 
attempt to effectively overcome this mishap by applying an 


easy trick. 
i:=0..rows(new_datamatrix02) - 1 
r* := cols(new_datamatrix02 ; ) 
corr_norm„diffj := 

T 

corr„norm_diff = 

This last matrix named "corr_norm_diff" is the one we 
intended to find. In each row of unequal element number, it 
contains a different random variable as function of our 


submatrix(norm_diffj, 0, rows (norm - 1,0,- l) 



0 

•1 

2 

3 

4 

5 

6 

0 

[1,32] 

[1,46] 

[1,51] 

[1,47] 

[1,50] 

[1,48] 

[1,48] 
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"devised” groups. Let us take a closer look in one of them, 
say the 5th. 

/ T \ <o> 

testj) := \corr„norm_diff ; 
testj) 0 = 

mean( testj) 0 ) = 0.197 
testj := 

testj o = 

mean(te$tJ 0 ) = 1.387 

THESE ARE % DIFFERENCE 

Defining the number of the bins: 

^points := 200 

ii := 0.. NpQjjjjs — 1 

j := 0.. rows(corr_norm_diff) - I 

upperj := ma>(corr_norm_diff j) 

lowerj := tnm(corr_norm_diff j) 

upper_max := ma>< upper) 

lowerjrnm:= min( lower) 

upperjnax - lowerjtum 

step := - 

N • 

iY pomt$ 

upperjnm = 14.215 
Iower_min = -10.339 
step = 0.123 

genjntervals „ := lower„min *f step ii 



rows (gen Jntervals ) = 200 


buildL(L) := L new *-L 0 T 

for q€ 1.. rows(L) - 1 
A 

^ stack ^k-,£ W , a) 


I^-buildUL) 



cols^) = 199 
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qq-O-Npdft-1 

avg_pdf qq := mean(L pdf ^ ) 


^avg_pdf = 0.998 


sen intervals = 



0 

1 

2 

3 

4 

"M 

6 

7 

0 

-10.339 

-10.216 

-10.094 

-9.971 

-9.848 

-9.725 

-9.603 

-9.48 


209 























































































































































































avg_pdf 



gjfl;. 

82 

83 .... 



ipiiilil 


0 

0.463 

5.015-10 -3 

6.468-10 *3 

4.614-10 -3; 

3.483-10 “3 

7.664-10 -3 

5.727*10 ”3 


avg_pdf g j := 0.005 
^Tavg_pdf =0.54 

avgjpdf 


avg_pdf := 


^avg_pdf 
^avgjpdf =1 


Npdfc-l 


mean_rv := ^ gen_intervals q o-avg_pdf 
qq=0 

step = 0.123 

m 


qq 


Npdsrl 


mean 


_sq_rv := ^ (genjntervals qq )"-avg_pdf 


qq 


qq=0 

mean_sq_rv = 24.219 

*9 

var_rv := mean_$q_rv - mean_rv~ 
varjrv = 15.658 
meanjrv = 2.926 

Defining Gaussian pdf: 


1, . ; . ' 2> var_rv 


. ' J; ; .;' feTsvaxjcv 
room^gain^ := genjntervals 
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END of MathCAD file used in calculating the excess 
gain/loss random variables (Third location) 

MathCAD file used in calculating the excess gain/loss 
random variables (Forth location) 

Calculating the normalized difference random variables 
random position STACK 

In this program we are aiming at calculating all the random 
variables of the normalized difference between the room's 
and the free space field strength (Etot and Etot_free 
respectively)-We start by reading the associated data files 
generated by NEC_3SC. 

NEC_B SC_dara 1 := READPRN'’C:/texnp/Thesis/x 1 c.txt") 

Ml := rows(NEC_BSC_datal) 

Ml = 213 

N1 := cols (NECJB SC_datal) 

N1 = 9 

NEC_BSC_data2:= READPRK”C:/temp/Thesis/x2c.txt’’) 

M2 := rows (NEC_B SC_data2) 

M2 = 133 

N2 := cols(NEC_BSC_data2) 

N2 = 9 

NECJBSC_data3 :=■ READPRK”C:/temp/Thesis/x3c.txt'') 
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M3 := rows (NEC_B SC_data3) 

M3 =49 

N3 := cols (NEC_BSC_data3>) 

N3 = 9 

NEC_BSC_data4:=READPRM;"C:/temp/Tfaesis/x4c.txt") 
M4 := rows(NEC_BSC_data4) 

M4 = 49 

N4 := cols ( NEC_B SC_data4) 

N4= 9 

NEC_BSC_data5 := READPRK"C:/temp/Thesis/x5c.txt") 
M5 :=rows(NEC_BSC_data5) 

M5 = 77 

N5 := cols (NEC_B SC_data5) 

N5 = 9 

NEC_BSC_data6 := READPRIsTC:/temp/Thesis/x6c.txt") 
M6 := rows(KEC_BSC_data6) 

M6 = 497 

N6 := cols(NEC_BSC_data6) 

N6 = 9 

NEC_BSC_data7 := READPRN"C:/tempyThesis/x7c.txt” ) 
M7 :=rows(NEC_BSC_data7) 

M7 = 133 

N7 := cols (NEC_BSC_data7) 

N7 = 9 

NEC_BSC_dataS := READPRiN["C:/temp/Thesis/xSc.txt") 
MS := rows(NEC_BSC_data?) 

MS = 49 

N8 := cols (NEC_BSC_data8) 

N8 = 9 

NEC_BSC_data9 := READPR-N( , 'C:/temp/Thesis/x9c.txt") 
M9 := rows (NEC_B SC_data9) 

M9 = 49 

N9 := cols(NEC_BSC_data9) 

N9 = 9 

NEC_BSC_datalO:= READPR^''C:/temp/Thesis/xlOc.txt") 
M10 := rows( NEC_BSC_datalQ 
M10 = 77 

N10 := cols (NEC_BSC_datal Q) 

N10= 9 

NEC_BSC_datal 1 := R£.4DPRN;"C :/temp/Thesis/x 1 lc.txt") 
Ml 1 := rows(NEC_BSC_datal I) 

Mil = 213 

N11 := cols (NEC_B SC_data 11) 

Nil = 9 
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3SIEC_BSC_data ^ 
ycoord ;= NEC_BSC_data (i) 
Zcoo^ := NEC_BSC_data ® 


P yi^cooii ^)” ~ (ycoord “ (^coord ~ l)” 

free_space_datal := KEADPRK"C:/temp/Thes;s/fs 1 c.txt" ) 
MM1 := rows(free_space_datal) 

MM1 =213 

NN1 := cols(free_space_datal) 

NN1 = 9 

free_space_data2 “READPRK"C:/temp/Thesis/fs2c.txt‘’) 
MM2 := rows(free_space_data2) 

MM2 = 133 

NN2 := cols(free_space_data2) 

NN2= 9 

free_space_data3 := READPKN("C:/temp/Thesis/fs3c.txt" ) 
MM3 := rows(free_space_data3) 

MM3 =49 

NN3 := co!s(free_space_data3 ) 

NN3 = 9 

free_space_data4 := READPRN;”C:/temp/Thesis/fs4c.txt" ) 
MM4 := rows (free_space_data4) 

MM4 = 49 

NN4 := cols (free_space_data4 ) 

NN4 = 9 

free_space_data5 := READPRI^"C:/temp/Thesis/fs5c.m") 
MM5 := rows(iree_space_data5 ) 

MM5 =77 

NN5 := cols(free_space_data5) 

NN5 = 9 

free_space_data6 := READPKK”C:/temp/Thesis/fs6c.txt” ) 
MM6 := rows(free_space_data6) 

MM6 = 497 

NN6:= cols(free_space_data6) 

NN6= 9 

free_space_data7 := KEADPRN[”C'7temp/Thesis/fs7c.txt") 
MM7 := rows(free_space_data7) 

MM7 = 133 

NN7 := cols(free_space_data7 ) 

NN7 = 9 

free_space_data8 := RE.4JDPRK"CVtemp/Theas/fs8c.txt" ) 
MM8 := rows(free_space_data8) 
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MM8 = 49 

NN8 := cols (free_$pace_data8 ) 

NN8 = 9 

free_space_data9 :=READPRH"C:/temp/Thesi$/fs9c.txt'' ) 
MM9 := rows(free_space_data9) 

MM9 = 49 

NN9 := cols(free_space_data9) 

NN9 = 9 

free_space_datalO := READPRN "C:/temp/Thesis/fs 1 Oc.txt" ) 
MM10 := row$(free_$pace_datalO) 

MM10 = 77 

NN10:= cols(free_space_datalO) 

NN10= 9 

free_space_datal 1 READPRK"C:/temp/Thesis/fsl lc.txt” ) 
MM 11 := rows(free_space_datal 1) 

MM11 = 213 

NN11 := cols(free_space_datal I) 

NN11 = 9 


( NEC_BSC_data ® ^ 




10 


20 


NEC_BSC_data 


<5> 


Ey := 10 


20 


NEC_BSC_data 


< 7 > 


Ez:= 10 


20 


Etot :=^) 2 -(Ev ) 2 + (E,) 2 
rows(Et o; ) = 1.539x 10 3 


/ < 3 > \ 

frce„spaoe_data_ 1 


^„free :=: 


10 


20 


free_space_data 


< 5 > 


Eyjtee- 10 


20 


frce_space„data 


< 7 > 


;_free’ : 


10 


20 


Etot_frec :— \/(Ex_fre^~ ” (Ey_fret) + (Ez_fre^~ 

Since the NEC-BSC dynamic range is only 30dB, we will have 
to account for this fact by "filtering" out all values 
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which are out of that range, 
fashion: 


Phis is done in the following 


Etot_dB (20-log(E tot )) 


ij:=0.. rows(E K)t _ dB ) - 1 

lowest_alIowed_Etot_value := - dynjrange 

lowest_allowedJStot_value = 28.39 

%ot_free_dB * = (lO-log^Eotj&e^)) 

lowe$t_allowed_Etot_fr_vaI := - dyn_range 

lowest_aIlowed_Etot_fr_vaI = 28.51 

Etotai dB.. :=ifjp^ 0 t_dB..j < lowest_allowed_Etor_value ,lowest_allowed_Etot_value 

E olal free dB := ifife* free dg < lowestjillowed_EtotJf_val, lowest_allowed_Etot_fr_val free 

~ “ jj \ ~ jj 2 

Converting to linear: 


^routI„dS 


^totai ^ 


20 


EmalfreedB 


20 


%otal_£ree-“ 

Let us now create the two datamatrices, namely the 
"datamatrix" which consists of the room's field strength 
values for every observed point and the "freespacedataOl" 
which includes the corresponding data for the free space 
case. 

datamatrix := augment(p , E<y ? i ) 
teespacedataOl := augment(p. E^i 

At this point, owing to the nonuniform layout of the 
distance (from the antenna) points in the □ array, we will 
have to distinguish them into several different groups 
using the following approach: 
lower := floor{min(p)) 
upper := ceil(ma^p)) 

waviTmrm and minimum distance for this antenna locations 

upper = 15 
lower = 0 

Groups separation for this calculation (in cm): 

separation :=50 

"Step factor" for this calculation: 


factor := ceil 


./ 100 


separation 

factor = 2 

N umb er of distance bins for this calculations 
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bins := upper-factor 

bins = 30 

upper - lower 

h :=- 

bins 

h = 0.5 


n := 0.. bins 


k n := lower + n-h 

It should be noted that the number of bins may be chosen in 
a fashion which best suits our desirable "distance steps". 
For instance, should we choose to have groups points 
separated by 50cm, where the maximum distance from the 
antenna is 49m, our intuitive selection should be 98. 

We next introduce a new function which shall be utilized 
several times in our subsequent calculations. As input, it 
expects a two column matrix and an array. It looks for 
those values of the first column (i.e., in our problem, the 
distance array □..which fall within our made-up groups and 
then creates a new matrix that contains all those second 
column's values whose respective first column value satisfy 
our requirement. It finally returns this new matrix which, 
it should be underlined, normally has a different number of 
elements in each row. 


build(M, bi nvector ) 


AeM^ 
B 4- M ^ 


k 4- binvector 


for le O..rows(k) - 2 
D; 4— 0 

for ie O..rows(k) - 2 
for je 0..rows(M) - 1 
D ; <— augment (Bj,Dj) if lq < Aj < lq +1 
for i € 0.. rows(D) - 1 

Dj e- submatrix(D i ,0,rows(D i ) - 1,0, cols (D;) - 2) if D; ^ 0 
D 

It is high time we used this helping function. 
datamatrix02:= build(datamatrix, k) 


ffeespacedataQ2 := build(freespacedata01 ,k) 


datamatrixO^ 




1 

2 

3 

4 

■}$()' 



[1,14] 

[1,30] 

[1,33] 

[1,63] 

[1,71]l 

[1,76] 

[1,94] 


The following new function is, seemingly, of high 
importance since it is the one which calculates the 
"normalized difference" random variables. It receives two 
matrices of the same dimensions and returns a new matrix 
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which in each row contains the values of the random 
variable pertinent to that row's "distance group." 
difference(X,Y) := 


for ke 0..rows(X) - 1 
C k <- 0 

(for js 0..rows(Xi - 1 
n <- cols(Xj) 

A 4— Xj • 

B 4— Yj 
i 4—0 

for i e 0.. n - 1 
dxffj 4- 20-logj 


A 0. i 

Bo. 




n 4- 


Cj ^—stacfc 


c j <_ 
k 4- rows(Cj) 

1 4— COis(Cj) 

Gj <- submatrix(Cj, 0, k - 1,1,1 - l) 


) 


The "room gain" is defined here in dB 

We now get rid of the redundant first row zeros and then 
calculate the normalized difference variables. 
new_datamatrix02:= submatrixfciatamatrixOZ i,rows(datamatrix02> - 2,0,cols(datamatrix.02 - 1) 
new_ffeespacedata02 :=$ubmatrix(freespacedara02,l,rows(ffeespacedata02) -2,0,cols(ffeespacedata02) - 1) 


norm_diff := difference(new_datamatrix02, new_freespacedata02) 
norm_diff r = 


Surprisingly and unexpectedly the Mathcad appears to suffer 
form a bug for it returns an erroneous matrix. The cause of 
this weird occurrence remains unknown, however we shall 
attempt to effectively overcome this mishap by applying an 
easy trick. 


n 


1 

IS 

« i *0 A 

. "4 

|i' .5 ” 


0 

[1,30] 

[1,33] 

[1,63] 

[1,71] 

[1,76] 

[1,94] 

[1,99] 


i := 0.. row$(new_datamatrix02) - 1 
r; := cols (new_datamatrix02 : ) 
corr_norm_diffj 

T 

corr_norm_diff = 

This last matrix named "corr_norm_diff" is the one we 
intended to find. In each row of unequal element number, it 


submatrix(norm_diff;, 0, rows (norm_diff;) - 1,0, r, - l) 



0 

'I 




izcm 

6 . 

0 

[1,30] 

[1,33] 

[1,63] 

f[1,71] 

[1,76] 

[1,94] 

[1,99] 
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contains a different random variable as function of our 
"devised” groups. Let us take a closer look in one of them, 
say the 5th. 

te$t_0 := 

tesiJ3 0 = 


mean(test_O 0 ) = 0.529 

/ j\ 0) 

test_J := \corr_norm„diff j 



0 

; i: V; 

2 vn 


■j//' ; 4';;; 


■44 : : 

7; ;■ 

8 

0 

0.171 

-4.62 

0.389 

-0.842 

-2.145 

3.178 

2.166 

-3.219 

1.823 

i_ 


mean(test_l 0 ) = -0-346 

THESE ARE % DIFFERENCE 

Defining the mnriber of the bins: 

•^points 

ii := 0.. NpQ inls — 1 

j := 0.. rows(corr_norm_diff) - 1 

upper j := ma>(corr_norm_diff j) 

lowerj :=mii^corr„norm_diffj) 

upper_max := ma>( upper) 

lower_min := min(lower) 

upper max - lowerjnin 
step := .. . 

1ST . 


( T \ ® 
\corr_norm_diff / 



0 

1 

2 

3 4 

;;;44 

6 

7 

0 

0.078 

-2.592 

0.235 

-0.191 

3.503 

3.511 

-10.335 

-0.102 


upperjnax = 12.092 
lower_min= -10.335 
step =0.112 

genjntervaIs „ := lowerjrnin + step *ii 



^ ; cols(coiT_nor m_ diffj) 


rows (gen..intervals ) = 200 


buildKL) := 


Lnew ^0 


for q € 1.. rows(L) - 1 
A«-Lq T 

Lnew stack(L ttew . A) 


Lpjf := buildL(L) 



cols(Lpdf) = 199 
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Antenna at location 4 













room^gain^B.. 


n 


END of MathCAD file used in calculating the excess 
gain/loss random variables (Forth location) 


End of all MathCAD files 
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